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L'integrale de chemin : 

L'integrale de chemin est une methode generale qui permet d'exprimer les prob- 
ability de transition en mecanique quantique, et plus generalement les fonctions 
de correlation en mecanique statistique et les propagateurs en theorie quantique des 
champs, a l'aide d'une somme sur les chemins possibles du systeme. Suivant le point 
de vue que Ton adopte, on peut y voir : 

• Une nouvelle formulation de la mecanique quantique, qui eclaire le rapport 
avec la mecanique classique et qui s'est revelee un outil tres utile, voire indis- 
pensable, en theorie quantique des champs, notamment pour les theories de 
jauge. 

• Une methode qui permet de formuler differents types de champ moyen, et qui 
fournit une approche systematique pour les developpements perturbatifs. 

• Un concept unificateur qui rend manifeste l'analogie tres etroite entre les 
phenomenes critiques et la theorie quantique des champs. 

L'objet de ce cours est de proposer une introduction a ces differentes applications 
de l'integrale de chemin. 
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Chapitre 1 



L 'integrate de chemin en mecanique 
quantique 

1.1 Introduction 

En 1948, Richard Feynman a initie une veritable revolution - meme si elle n'a 
pas ete soudaine et est encore en marche - en proposant une nouvelle formulation 
de la mecanique quantique. L'idee de depart, due a Dirac, est la recherche d'une 
formulation de la mecanique quantique directement en termes du lagrangien. La 
formulation hamiltonienne n'est en effet pas completement satisfaisante pour au 
moins deux raisons : 

dL 

• Elle suppose l'existence d'un moment conjugue Pi = — — pour chaque variable, 

% 

ce qui n'est pas le cas le plus general (c/les difficultes de cette methode pour 
le champ electromagnetique). 

• Elle n'est pas manifestement covariante par une transformation de Lorentz. 

Pour mieux comprendre les motivations de Dirac et Feynman, il est utile de faire 
un bref rappel des differentes formulations de la mecanique quantique. 

1.1.1 La theorie des quanta (Bohr-Sommerfeld, 1915) 

La premiere theorie de la quantification des niveaux d'energie dans les atomes due a 
Bohr et Sommerfeld (1915) est basee de fagon tout a fait essentielle sur la formulation 
hamiltonienne de la mecanique classique. Les etapes sont les suivantes: 
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Determiner le lagrangien L(q 1 , . . . , q N , qi, . . . , q N ) (suppose independant du 
temps). 

En deduire le hamiltonien par une transformation de Legendre: 

H(q 1 ,...,q n ,p 1 ,...,p N ) = J^Pifc- L{q 1 ,...,q n ,q 1 ,...,q N ), (1.1) 

i 

dL 

ou les Qj sont des fonctions de (qi, . . . , q n ,Pi, ■ ■ ■ ,Pn) donnees par = -—. 

' ' % 

Resoudre l'equation de Hamilton-Jacobi 

et en deduire une transformation canonique vers les variables action-angle 1^ 
cui, telles que: 

( dl 

dt 

(1.3) 

dwi 

I 

Les regies de quantification sont alors donnees par: U = hni, ni — 0,1, ... . Pour un 
systeme a un degre de liberte, cela s'ecrit: 

/ = — <j> pdq = nh, (1-4) 

h 

ou h = — (h— constante de Planck) = quantum d'action. 

Cette theorie, consideree par Sommerfeld comme "la voie royale vers la quantifi- 
cation", s'est en fait rapidement revelee insuffisante. Sa limitation majeure est de 
reposer sur une solution complete de l'equation de Hamilton-Jacobi. Cela suppose 
que le systeme est separable. Or les systemes separables constituent l'exception: 
l'immense majorite des systemes mecaniques ne sont pas separables, et il s'est avere 
impossible de generaliser cette formulation sans introduire de nouveaux concepts. 



1.1.2 La theorie des matrices (Heisenberg, 1926) 

En 1926, Heisenberg a demontre que les memes regies de quantification pouvaient 
etre obtenues pour un systeme separable si les energies etaient definies comme les 
valeurs propres de l'operateur hamiltonien H: 

v 2 

H = — + V(x) (1 particule), (1.5) 

2m 
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ou x et p doivent etre consideres comme des operateurs satisfaisant les regies de 
commutation 

[x,p} = m (1.6) 

Mais cette formulation, qui peut etre etendue automatiquement a n'importe quel 
systeme comportant plusieurs particules en introduisant le hamiltonien 

2 

H = J2^ + V(x 1 ,...,x N ) (1.7) 



et les regies de commutation 



\pi,Pj] = (i-8) 



ne necessite pas que le systeme soit integrable. 

Heisenberg l'a done proposee comme une theorie generale de la quantification. C'est 
effectivement la premiere theorie de la mecanique quantique. L'idee d'introduire 
des operateurs qui ne commutent pas etait basee sur les proprietes des matrices, 
et Heisenberg a utilise cette representation pour determiner les valeurs propres du 
hamiltonien. 

Cette fagon de quantifier le mouvement, connue sous le nom de "quantification canon- 
ique", repose toujours de fagon essentielle sur la formulation hamiltonienne du prob- 
leme de mecanique classique correspondant. En particulier, revolution temporelle 
des operateurs est regie par l'equation 



dt 



0,H 



(1.9) 



1.1.3 La mecanique ondulatoire (Schrodinger, 1927) 

Independamment de Heisenberg, et en se basant sur les proprietes d'onde des elec- 
trons revelees dans les experiences de diffraction, Schrodinger a formule une autre 
approche basee sur le concept de fonction d'onde. Dans cette approche, l'etat d'un 
systeme est repere par une fonction d'onde dont revolution temporelle est regie par 
l'equation 

d k 2 d 2 

i^-^(x, t) = - ^ a^(*> *) + V(x)if>(x, t) (1.10) 

(pour 1 particule non relativiste). 

Les niveaux d'energie d'un atome correspondent aux etats stationnaires 

ih^(x,t) = E*(x,t). (1.11) 
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Dans cette formulation, le hamiltonien n'apparait pas explicitement. Cependant il 
est possible de faire le lien avec la formulation de Heisenberg en remarquant que 
l'equation peut se recrire 

ik^ = Hip, (1.12) 

p 2 d 

avec H = h V(x) et p = —ih—. 

2m ox 

En effet, les operateurs x et p definis par 

x : I— > x(f> 

rlrh 

p (f) i — * —ih 



dx 



satisfont [x,p\ = ih. 
Demonstration: 



[x,p}(j) = x(-ih)^ + ih-^(x(f)) 

= x(—in)— — \-ihx— — h lhtp 
ox ox 

= ih(f). 



Le lien entre les deux representations en termes de bases de l'espace de Hilbert a 
ete etabli par Dirac. 



1.2 L'integrale de chemin (Feynman, 1948) 



L'evolution temporelle de la fonction d'onde peut se decrire a l'aide de l'operateur 
d'evolution: 

m = u(t,t )^(t ) (i.i3) 

ou U(t,t ) est un operateur qui satisfait l'equation 

dU 

ih^- = HU. (1.14) 
dt v ' 

L'operateur d'evolution peut etre aussi defini par ses elements de matrice dans une 
base de l'espace de Hilbert des fonctions d'onde. Deux bases sont particulierement 
utiles: \x) et \p). 
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\p ) Les kets \p ) sont definis comme les vecteurs propres de l'operateur p de valeurs 

' d 
propres p. Or, p = —ih—. La fonction d'onde (j) p (x) associee a p satisfait done: 

d 

Cette equation est manifestement satisfaite par e ipa:7 \ La seule subtilite reside 
dans la normalisation. Comme e ipx ^ n n'est pas de carre sommable, on ne peut 
pas la normaliser a 1. La convention de normalisation est fixee par la relation 
de fermeture: 

J dp|p) (p| = l (1.16) 



Posons 



\P) 



_L e iWM (L17) 



J dp\p) (p\(f)(x)) = ^ J dp J dx'e ipx/h e- ipx ' /h (f)(x'), 

= J dpe ipx/h J dx'e- ipx ' /h <j)(x'). (1.18) 

Or, la transformee de Fourrier de <f>(x) est definie par 

4>(k) = -±= [ dxe^^x), (1.19) 

et la transformation inverse est donnee par: 

4>{x) = -jL [ dke ihx 4)(k). (1.20) 
ylix J 



J dp\p){p\4>{x)) = j^J dpe^/»^(| 

= ^/ ^ x m = ^(x) (i.2i) 

N = (1.22) 



\x) Les kets \x) sont les vecteurs propres de x. Or, par definition de l'operateur 
x, on a: 

x(p(x) = xcj)(x). (1.23) 

Une fonction propre <p Xo de l'operateur x de valeur propre x Q satisfait par 
ailleurs: 

x(fr X0 (x) = x <fr X0 (x). (1.24) 
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<p Xo doit done satisfaire l'equation: 

x(p XQ (x) = x (p xo (x) pour tout x. (1-25) 

On en deduit que <f> XQ {x) — si x ^ xq et que 4>xq(%o) psut prendre n'iniporte 
quelle valeur. Mais pour que cette fonction n'ait pas une integrate nulle, il faut 
qu'elle soit proportionnelle a la "fonction 5" de Dirac 5(x — x ) definie par: 

J 5(x — Xo)f(x)dx = f(x ) pour toute fonction /. (1-26) 

Pour la normalisation, on va encore imposer une relation de fermeture. Or, le 
choix 

\x ) = \8(x-x )) (1.27) 
conduit a la relation de fermeture 

dx\x)(x\ = l. (1.28) 

En effet, (<p\x) = (jf(x) et (x\<j>) = (j>(x). Ainsi, 

(4\x)(x\4T) = <P{x)<V{x) 
Jdx{<f>\x) {x\(f)') = J dx(f)*(x)(f)'(x) = 

C.Q.F.D. 

Par ailleurs, la transformee de Fourier de la fonction S(x — x Q ) est une onde 
plane: 

f 8{x - x ) e~ ipx/h dx = e - ipxo/h 



I 



La transformation inverse conduit a la relation suivante: 



J e ipx/n e -ipx /h dp = ^5{x - x ). 



Le propagateur (ou fonction de Green) est defini par 

G(x,t;x ,t ) = (x\U(t,t )\x ). (1.29) 

Si Ton connait le propagateur pour tout x et tout xo, on connait l'operateur d'evolu- 
tion, i.e. la solution de l'equation de Schrodinger. 

L'objectif de Feynman etait de relier le propagateur au lagrangien du systeme clas- 
sique correspondant. Or, la solution de l'equation de mouvement pour l'operateur 
d'evolution s'ecrit: 

U(t,t )=e-^ H , (1.30) 
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une expression qui fait intervenir le hamiltonien et non pas le lagrangien. 

L'idee de base vient de la remarque suivante: rhamiltonien H est la somme de deux 
termes H = T + V. 

• V(x) a pour etats propres \x ) : V(x)\x) — V(x)\x) 

p 2 p 2 

• T = — a pour etats propres \p) : T\p) = \p). 

2m 2m 

Du coup, une valeur moyenne du type ( x|e~ AT e~ Al/ |x ) est facile a calculer. II suffit 

d'inserer la relation de fermeture 11 = J dp \p) {p\ (on a pose 

- WWW - Jm^\pHp\^M 

= J dp(x\p)e- x ^(p\x )e- xv(xo) . (1.31) 
Malheureusement, comme T et V ne commutent pas, 

e -HT + V)^ e -XT e -XV_ (132) 

Comment faire? 

L'idee de Feynman consiste a utiliser ce qu'on appelle la formule de Trotter qui 
stipule que 

(x\e- x ^\x )= lim (xle-^e-^.-.e-^e-^lxo). (1.33) 
La demonstration procede en 3 etapes: 



• o" A ( T+y ) = ( e -* T+v V N ) N (evident) 

. e -X(T+V)/N = e -\T/N e -XV/N + q 



1 

N~ 2 J' 



f 2 

Eneffet, e e(A+B) = 1 + e(A + B) + -{A + B) 2 + . . . 



et e eA e eB = ( 11 + tA + ^-A 2 + (l + eB + ^-B 2 + ... 



2 

A 2 

2 

= 1 + eA + eB + 0(e 2 ) 

e e(A+B) = e eA e eB + ^ 



-XT/N e -XV/N^ N _ ^ e -X(T+V)/N^ N 



= o 
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Pour cela, on remarque d'abord que 1'identite 



x 



N -y N = (x-y) (x N ^ + x N ^y + ■■■ + y N ^) 



pour x, y G IR doit etre rearrangee pour deux operateurs qui ne commutent 
pas: 

x N -y N = x N '\x -y) + x N ~\x -y)y + -.- + (x- y)y N ~\ 

Dans le membre de droite de cette egalite, tous les termes sont de la forme 
x p y q , et ils s'annulent 2 a 2, comme dans l'identite habituelle, sauf x N et y N . 

Ainsi, {e~ XT/N e~ xv/N ) N - ( e - x{T+v)/N ) N est la somme de N termes qui con- 
tiennent tous e - XT / N e -* v / N _ e ~MT+v)/N _ Comme ce f acteur est 0(1/N 2 ), la 
somme est d'ordre 1/N. C.Q.F.D. 

Revenons a l'expression du propagateur: 

G(x,t;x ,t ) = lim ( x le-^e"^ . . . e'^^e"^ I x ) 

et du type J dx\x){x\ entre e XV ^ N et e XT / N . \\ vient: 

lim / dp N e~^N(x\p N ) / da^-i e" A1/(:Ejv - l)//7V (pAr|:rAr_i ) 
x J dp N - 1 e- X -^r{x N ^ 1 \ PN _ 1 ) J dx^e-^^-^ipM^lxN^) 



G(x,t;x ,t ) = 


lim / 

Af^+oo J 


X 


J &Pn-i 


X 




X 


y dp 1 e" 



„2 



^(^ibi)e- Ay(:ro)/Ar (pi|xo) (1.34) 



L'expression devient done l'integrale du produit de facteurs de la forme 

J dp(x\e- XT/N \p) (p\e~ xv/N \y) = J dpe'^ (x\p) e' xv{y)/N (p\y). (1.35) 



Mais (x\p) = 1 e iW^ L'integrale devient: 

e -XV(y)/N f 



2nh 

C'est une integrate gaussienne (voir Appendice A). 



J dpe-^e* p(x - y) . (1.36) 
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i I / 27T J 2 




Or, z(j)= dxe * Ax +jx = \ -re^ 



e -XV(y)/N , e -XV(y)/N l 2nm N (*-y)* m N 

2nh J F 2nh V A 

Cette expression reste valable si A est imaginaire pur a condition d' interpreter \fi 
comme exp(i7r/4). On en deduit l'expression suivante pour G: 



N/2 



f ( mN V 

G(x,t;x ,t ) = lim / dxi . . .dx^-i - _ . 2 

m(xj+i — Xj) 2 N ^V(xj 



N-l 



x | J exp 

j=0 



2Xh 2 



N 



avec xn = x. 
(1.37) 



(f-t ) ^A , iNe n 
Posons e = — — — = — => A = — — . On a: 
N iN h 



N 



Nh 



Xh 2 h 2 \Ne he 

A ie 

~N = ~h' 

Le propagateur peut done etre reecrit: 

// m \ N / 2 
dxi . . . dx N _i 
\2-K\heJ 



(1.38) 



N-l 



N-l 

L'expression e 

j=0 

it - 1 ) 



X 



Tft I Xj+i — Xj 

~2 



CXP ^ 



3=0 



V( Xj ) 



- \/(x, 



(1.39) 



rappelle une somme de Riemann. En 



effet, e 



N 



Considerons une fonction derivable x(s) definie entre t et i, et 



exprimons l'integrale 



ds 



to 



= (S'-VW.)) 



sous forme d'une somme de Riemann: 

ft 



/ ds 

Jt 



N-l 



in 



x 2 - V{x) 



= lim — — — — Xa — V(xj). 

j=0 



(1.40) 
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Mais Xi 



(t-t )/N' 



[ ds 

J to 



rn 



x 2 - V(x) 



iV^+oc \ N J ^— ' 



m / Xj+i — Xj \ 2 



V{x 3 



2 \(t-t )/Nj ' X ~ J ' 

(1.41) 

/ ds L(x, x, s) est Taction 
le long de la trajectoire x(s). Ceci suggere d'ecrire le propagateur sous la forme: 



m 



Or, '-^-x 2 — V(x) n'est rien d'autre que le lagrangien, et 




ou le symbole / Vx designe la somme sur tous les chemins possibles, modulo une 
mesure qui reste a definir. Le propagateur apparait done comme la somme de 
exp(j:S[x]) sur tous les chemins possibles du systeme (voir Figure 1). Cet objet 
s'appelle une integrate de chemin ou integrate fonctionnelle. Cette ecriture est tres 
suggestive, mais elle est purement formelle, et pour lui donner un sens, autrement 
dit pour definir la mesure f Vx, il faut imperativement revenir a la forme explicite. 
C'est d'autant plus vrai que le chemin typique n'est pas differentiable: comme Xj et 



Xj+i prennent toutes les valeurs entre — oo et +oo, 



Xj+i 



ne tend pas vers une 



constante mais diverge quand iV — > +oo (voir le mouvement brownien). 



x„=x 




Figure 1: Exemples de trajectoires qui interviennent dans la somme sur les chemins. 



Ceci etant dit, cette formulation repond a la question que s'etait posee Feynman (et 
avant lui Dirac) sur la possibilite de formuler la mecanique quantique a l'aide du 
lagrangien. Avec la definition de l'integrale fonctionnelle donnee par la correspon- 
dance entre les 2 expressions precedentes du propagateur, la mecanique quantique 
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peut etre basee sur cette definition du propagateur. En d'autres termes, on peut 
faire l'economie de la quantification canonique et prendre cette expression du prop- 
agateur pour quantifier le mouvement d'une particule dans un potentiel. 

La difficulty essentielle de cette methode, c'est qu'elle ne fournit pas de fagon interne 
la definition de l'integrale fonctionnelle. II faut par exemple rajouter une prescription 
pour obtenir une expression qui redonne l'equation de Schrodinger en presence d'un 
champ magnetique. De raeme, en theorie quantique des champs, le sens qu'il faut 
donner a l'integrale se refere a la quantification canonique. 

Notons que la generalisation a N degres de liberte est sans probleme: les integrales 
intermediaries sont des integrales de volume a N dimensions: dx { — > d N Xi. 

Exemple: Particule libre (V — 0) 

Dans ce cas, le propagateur se calcule immediatement: 

( x\e~ XT \xo ) = J dp (x\e~ XT \p) (p\x ) 

Ap 2 e ipx/h e -ipx /h 

dpe~ 

dpe~^ p2 e ipix ~ xo)/h 



= — I 

2nh J 



X . i(x-x ) 

A = — > 3 = * 

m h 

(x — xo) 2 m 



1 /27rm 



h 2 2X 



soit 



x 



i(t-«o) f, 2 

h 2m 



X 



m 



2nhi(t — t c 



cxp 



lm 



2k(t - 1 ) 



(x - x f 



Or, pour une particule libre, Taction se calcule aisement. La trajectoire classique 
est donnee par: 

x(t) =x + ^^(x-x ) (1.42) 

t — Iq 



S c \ = I drL(x,x) = I dr]-m( 
Jt Jt 2 \ 



2 \drj 



k 2 \t-to) 



1 (x- x ) 
-m 

2 t-t 



m 



lSr. 
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1.3 Oscillateur harmonique 



Considerons le cas ou V(x) = ^muj 2 x 2 . Le propagateur s'ecrit: 



G(x,t;x ,t ) 



lim ( - ., | 2 / dxi...dx 



N 

ie x - 



3=0 



Tfi I %jj r \ Xj 

~2 



2 2 

-moo X; 
2 J 



avec xn + \ = x et e 



t-t 



(1.43) 



_r e£ape: il faut travailler sur l'exponentielle. 

Tout d'abord, il faut isoler les termes qui contiennent x et x^+i puisqu'on n'integre 
pas sur xq et xjv+i- 

(Xi - Xq) 2 = x\+ X 2 Q - 2XqXi 

(x 2 — X\) 2 = x\ + x\ — 2xix 2 



{x N+ i - a; at) 2 = x 2 N+1 +x 2 N - 2x N x N+1 
L'argument de l'exponentielle peut done s'ecrire 

(2xj H 1" 2rc Ar _ 2x1X2 2x N _ l x N ) - 1 ^mw 2 (xl H h x^) 



+^-(-2^1) + ^-(-2x^x^+1) + ^7^0 + ^7v+i) - T~> ,h 
Cela se met sous la forme 



1m 



1m 



iel 



2^2 



( I N N \ im 

exp (^~2 + ex P ^ + 



avec 




ie 



/1 



+ — moo 

h 







Q 1 im 

et j = - — 



1 2/ 



/ a; 




V XN+1 j 
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ou on a laisse tomber le facteur exp 
parait a la limite N — > +00 (e — > 0). 

Or, 



le moo 2 



qui vient du potentiel car il dis- 



/ x JV N \ 

dxi... dx N exp I -- ^ XiAijXj + ^ jjX; j 



(voir formulaire). 



Comme seuls ji et sont differents de 0, on a 



(^4 



(A-XJi + iA- 1 )^* 



Posons 



A N = 




et calculons (A , (A 1 ) 1N , (A 1 ) N1 , (A 1 ) NN en fonction de det A N . 
De fagon generale, on a: 



7 j* 



« det A 



ou a u = cofacteur = (— 1)* +J det(-By), ou .Bjj est la matrice qui s'obtient a partir de 
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A en supprimant la ligne i et la colonne j. 



i 11 = det 



IN 




(-1)^+1 det 



det(A N _i) = a 



NN 



N+l iJV-1 . Nl 



(-l) ly + L b 



a: 



N-l 



^ A ~^ = TT^ - W det (^-i) +j%det(A N . 1 ) + 2 3l3N (-l) N+1 b N - 1 ) 
det Aw 



Par ailleurs, on peut aisement calculer det(A N ). En effet, on a 



det(A N ) = adet(A N _i) — bdet 



( b b — 

a b 



\ 





Ob a b — 



V 



adet(A N ^ 1 )) - b 2 det(A N ^ 2 )- 



Posons In = det A N . On a done une recurrence double qu'on peut mettre sous 
forme matricielle: 

In \ = ( a -b 2 \( I N -i 
In-i U Oil I N _ 2 



soit 



In 
In-i 



a -b 2 
1 



N 



Io 

1-1 



I et I -i sont fixes par la condition Ji = a, I 2 = a 2 — b 2 
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Pour calculer 



a -b 2 
1 



N 



= C N ', on diagonalise C 
A , 



c = s 



A_ 



s-\ 



A i et A_ sont solutions de 



det 



a - A -b 2 
1 -A 

A 2 - aX + 6 2 








A 



a ± Va 2 - Ab 2 



± 



La matrice S peut se mettre sous la forme: 



A+ A_ 
1 1 



Ainsi, 



A+ 








A_ 


A+ 








A_ 



En effet, det S 

et S- 1 



s- 1 = 



A+ - A. 



A+-A_ V -1 A 



1 



A. 



-A+A_ 



s- 1 



A+ - A_ V "A- A + A_ 

1 / A 2 + - A 2 -A 2 + A_ + A+A 2 
A + -A_V A +-A_ -A+A_ + A + A_ 

A + + A_ -A+A_ 
1 



Mais d'apres l'equation satisfaite par A+ et A_, on a: 

A + + A_ = a 
A+A_ = b 2 



C.Q.F.D 



On a done: 

In 
In -i 



S 



= s 



= s 



s 



-1 



A^ 

A^ 

A^ 

\ N ) X._._ A 

1 



A+- A_ 



-X N /d 



f \ \ 

y_ 
\ d j 
i 



A+ - A_ 



\ N + 1 — \ N + 1 
V - \ N 
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In = 



\ N + 1 — \ N + 1 
A+-A_ 



Revenons au propagateur. II s'ecrit: 



,. / m \—r- [ im , o 2 x 



x(2tt 



.AT/2. 



: eX P \ n 



1 — m 



2 ftV 



l + xl +l ) I ^+2x,x N+1 (-l) N+ ^ 
In 1 



N 



II faut desormais calculer cette limite. Les donnees sont: 



2m e 

a = I — — h -mto" I i 

he h 

he 



a ± Va 2 - Ab 2 t - t 

A ± = 7; > e 



N+l N+l 



1. Calcul de A + et A_: 



6 2 

a 2 - 46 2 



4m 2 m 2 w 4 e 2 4m 2 w 2 
+ 



m 



ft 2 e 2 
2 



ft 2 



h 2 



2.2 



re 

Am 2 oj 2 m 2 oj 4 e 2 



h 2 



A + = 



ft 2 

-m ?nw 
-1 + 



+ 0(e) 



he h 



. . 2muj 
K-\- = — + 0(e) 
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2. Calculde - \ N+l : 



\ N + +1 



m muj\ N + 1 



-i + 



he h 

■ \ N+l 

-mi 



•\ N+l 

-mi 



he 



1 + 



muo 



he \ 



N+l 



h —mi 



1 



\ N+1 



y^N+l _ yN+l 



he 
—mi 
~he~ 



—mi 



he 



1 + ieuj) 



N+l 



N+l 



N+l 



1 + 



1 - 



iuT 
N+l 

iuiT 
N + l 



N+l 



N+l 



■\ N+l 

-mi 



he 



1 + 



iooT 
N+l 



N+l 



N + l 



N+l 



3. Coefficient devant les exponentielles: 



\2mheJ v 1 



V27ri^e/ v ; 



^W+l _ \ N + 1 

n {2muj)^ 



X 



he 



-mi 



N+l 



1 + 



iooT 



N + l 



N+l 



- 1- 



N+l 



N+l 



. n+i.n ( 2mu\ 2 



N + l 



N+l 



/mcu\ 

\~^h) 



1 + 



iooT 
N + l 



N+l 



iuoT 
N + l 



N+l 



iooT 
N + l 

-1/2 



-1/2 



N+l 



-1/2 



Mais lim ( 1 + 

N^+oo V AT 



= e x . En effet, posons -ujy = (l + j^j 

In u N = Nln (l + — ) = N ( — - - 

N \ N) \N 2 



N 



+ ... 



X. 



Ainsi 



lim 



1 



iu)T 



N+l 



iooT 



N+l 



N + lj \ N + l, 

Le coefficient devant les exponentielles est done: 



-1/2 



= [2isincjT]- 1/2 . 



muj 



2ii\hsiYiojT 
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4. Coefficient de x xn +1 : 

1_m2 2( 1)^+1 6JV-1 ^ 1 ~ m \( ±)N+i fm.\K-i 2mu ( In 



2h 2 e 2 In 2 h 2 e 2 \he J h I -mi 7 2isinu;T 



mu i / 1 

-muj- 



ihsinuT h \ sin uT 



5. Coefficient de x 2 

N 



'0- 



In \ -mi 



fte \ 1 



h 4- jiHT_\ N+1 _ (l _ ™T 
\ l ^ N+l) \ L N+l) 



he \ sincoT(l- ±) 



—mi / sin wT 

Tie \ sin ouT cos(uT / N) — sin(uT/N) cosuT 



—mi J sin cjT 

Tie he cos u;T . _ , „ T , 
- + :- — sm(uT/N). 



—mi —mi sin cuT 

Utilisons cette expression pour calculer le coefficient total: %1 est facteur de 

im 1 m 2 In-i im 1 m 2 f he he cosuT . \ 

We~2'W?~l N ~ = ^e~2¥?\^i + mi^Jf Sm{ 1 ') 

im 1 mi 1 mi cos uiT , m . „ T , 

= t ^ + ^ sin(uT/N) . 

2he 2 he 2hesinuT K 1 J 



Le dernier terme tend vers 



i moo cos uT 
h 2 sin luT 



puisque 

i sin(cc>T/iV) = sin(cuT/N) -> w. 

Recapitulation: 



mw 



2sincuT 



((a^ + Xf) cos luT — 2xixf) 



Interpretation: Afin d'interpreter ce resultat, calculons Taction classique de l'oscillateur 
harmonique. 



1 . 2 mw 2 
-mx — . x 

2 2 



mi + mu x = 



# = #i cos o>T + X2 sin cuT 
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x 



f 

Jo 



x 2 dt = 



t = 



t = T 



x = x f 



X 2 



X\ Xi 

Xf — Xi coscuT 
sincuT 



Jo 



dtL(x, x) 



dt 



mcu 



-mx — 



-x 



x\ cos 2 cut + x\ sin 2 cut + 2x\X 2 cos cut sin cut 



2 cos 2cut +1 2 '. 

= 2^ - h^2" 

T 



cos 2cut 



sin 2u;t £ 
Acu + 2 



+ x 2 



o 



+ X\x 2 sin2cjt 

T 

- x x x 2 



sin 2cjt 



/ sin 2cuT 



2cu 



+ T) + 



T — 



Acu 
sin 2T 

2cu 



o 



cos 2out 

2cu 
cos 2o>T - 



- x ± x 2 - 



2cu 



J" 

Jo 



x = —Xicu sincut + x 2 cu coscut 



x 2 dt 



xicu 2 ( sin 2cuT \ x 2 cu 2 



2cu 



+ T\ + 



+X 1 X 2 0j' 



,cos2cuT- 1 



2cu 



T — 



sin 2cuT 
2cu 



Finalement, 



S 



m 



2 2 
a; x 



)dt 



m 
2~ 



sin 2u;T x?a; sin 2cjT 2 cos 2a; T - 1 

— + x x x 2 cu 

2 cu 2 cu cu 

\-x\ sin 2cuT + x\ sin 2cuT - Ax x x 2 sin 2 cuT] . 



Reexprimons le terme entre crochets [. . .] a l'aide de x,- L et xf. 

x 2 + xf cos 2 cuT — 2xix f cos cuT 



x\ 



Xi 

x±x 2 



sin 2 cuT 



Xi 

XiXf—x? cos u>T 
sin u)T 



[...} 



X.: 



2 sin ouT cos cjT + 
coscuT 



+4 sin 2 cuT 

+XiXf 



2 sin cuT cos cjT cos 2 cuT 
sin 2 cuT 
. sin cuT cos cuT 



sincuT 
2 cos cuT 
sin 2 cuT 



+ x1 



sin 2 wT 



2sincuT cost-XT -4 



sin 2 cuT 
sin wT 
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xj 



X, 



XiXf 



-2 sin ujT cos ujT + 2 



cos 3 ujT 
sin ojT 



+ 4 sin u/T cos u;T 



sincuT 
1 



[—2 cos ujT sin 2 cuT + 2 cos 2 cuT cos cjT + 4 sin 2 cjT cos ojT\ 



sin wT 

cos cuT 

= 2 - 

smcuT 

_ cos ujT 



[2 cos wT(sin 2 uT + cos 2 wT)] 



sincuT 
1 



sinwT 
-4 

sincuT 



[-4 cos 2 loT - 4 sin 2 cuT] 



Finalement, on trouve 



S, 



muj 



cl 



2sincjT 



[(a; 2 + x 2 ) cos cjT — 2xiXf] . 



Conclusion: 



cxp 



h 



Comme dans le cas d'une particule libre, le propagateur est etroitement relie au 
mouvement classique. La phase est egale ai/k multiplie par Taction classique. 



1.4 Approximation semi-classique 



Dans le cas d'un potentiel general, on peut faire les deux remarques suivantes: 



• L'integrale fonctionnelle ne peut en general pas etre calculee exactement. 

• Le resultat obtenu dans le cas de l'oscillateur harmonique suggere que le por- 
pagateur peut etre approche par 



G(x, t; x , t ) — cte x e 



k S cl 



(1.44) 



Est-ce effectivement une bonne approximation? 
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Reponse: si le chemin classique existant est un minimum bien defini de Taction 
classique (voir plus loin pour plus de details), alors c'est une bonne approximation. 
C'est en fait le terme dominant dans un developpement en puissances de h. 

L'idee de base est de generaliser la methode de la phase stationnaire a l'integrale 

ftf 

de chemin. En effet, Taction S[x] = / dtL(x,x) est extremale pour le chemin 

Ju 

classique, et une approximation du type phase stationnaire devrait conduire a une 
expression contenant e^ 5cl . 



Considerons done Taction 



S[x] 

avec L(x, x) 



/ dtL(x,x) 
Ju 



-mi 2 — V(x) 



5S[x] 
5x(ti} 
5 2 S[x] 
5x(ti)5x(t 2 ) 

=>■ S[x + rf\ 



dt [mx5\t - h) - V'(x)S(t - h)} 

dt [m5'(t - h)5'(t - t 2 ) - V"(x)5(t - t!)S(t - t 2 )} 



f f f , 5S[x] 
Ju 



Sxih)' 



\ j ^ J dt^Mt^ 



S 2 S[x] 



)5x(t 2 ) 



+ ... 



La trajectoire classique satisfait 

5S[x] 



Sx(ti 

En effet, peut etre explicite en 

5S[x 



= (Vti). 



X — X c \ 



5x{ti 



= -mx(h) - V'ixih)), 



(1.45) 



(1.46) 



et d'apres Tequation d'Euler-Lagrange, le membre de droite est nul. 

S[x d + r]] = S[x cl } + ^ J dt J dh J dt 2 vihMh) (1.47) 
x [ m S'(t - tJS'it - t 2 ) - V"(x cl )S(t - tJSit - t 2 )} (1.48) 



S[x cl \ + J dt 



^vW - \v"{x c Mtf 



(1.49) 
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Le deuxieme terme est Taction d'un oscillateur harmonique dont la frequence depend 

muj 2 .. ... 
du temps: — — = V \x c \{t)). 



Par analogie avec 1' approximation de phase stationnaire, l'etape suivante consiste a 
decomposer l'integrale de chemin de la fagon suivante: 

J Vxe* s[x] = e* s[x ° l] J T>rie^ s[xd+v] - s[x * ]) (1.50) 

ou dans le membre de droite la somme est limitee aux chemins rj tels que rj(ti) = 
r)(tf) = 0. Si Ton se limite au developpement au second ordre de S[x c \ +rj] — 5[x c i], 
on arrive a l'expression 

G(x, t; x , t ) ~ e^G(0, t; 0, t ) (1.51) 

ou G est le propateur d'un oscillateur harmonique dont la frequence depend du 
temps. 

Remarques: 

• Dans la mesure ou l'integrale de chemin n'est definie que par la procedure 

lim , manipuler Taction n'est pas legitime a priori. C'est toute la puissance 

mais en meme temps Tambiguite de l'integrale de chemin. Elle suggere des 
manipulations tres efficaces mais difficiles a justifier a priori. On essaie plutot 
de verifier la validite du resultat a la fin du calcul. 

• Pour Toscillateur harmonique, 

5 3 S[x] 



. — = et V" = mo? = cte (1.52) 

<kc(ti)&c(f 2 )£a;(*3) 

Le developpement s'arrete done au 2eme ordre, et on devrait retrouver le 
resultat exact. C'est effectivement le cas. Le calcul que nous faisons conduit 
a un propagateur egal a 

e H*ci] x G(0,T;0,0) 

puisque rji — rjf = 0. Or, ce propagateur est un cas particulier du propagateur 
deja calcule pour Toscillateur harmaonique, et il redonne le bon prefacteur. 
La methode est done "justifiee". 

Dans le cas general, le propagateur G(0,t;0,t ) est defini comme toujours par la 
limite d'une integrale multiple: 



G(0,t;0,t ) = lim [ dx ± . . .dx N (-^—) 



N+l 
2 




m 
2~e 



Xj+i - Xj) 2 - -ec 



3 X j 



(1.53) 
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avec Cj = V"(xj) 



Comme pour l'oscillateur harmonique standard, on peut ecrire l'exponentielle sous 
la forme 



( if A \ 

\ i,j=l J 



avec 




\ 



le 

+ n 



V c N J 



-i 

-1 2/ 

II n'y a pas de matrice j car rj(ti) = rf(tf) = 0. On en deduit 



G(0,t;0,t ) 



lim 



/ m \ — 
V2&J 



^ (2tt; 



AT 
2 



Vdet A 



( m \ 2 



lim (— — ) 



(^) det A 



l -i 

/ m \ 2 1 



lim I — ) 



VdetB 



avec 



5 





V 



(1.54) 



-1 2/ 

Posons Pn = det 5. Par le meme raisonnement que pour det A, pjy satisfait la 
relation de recurrence 



avec 



p N+1 = 2 Cat +1 J9at - 

m 



pi = 2 ci et p = 1- 

m 



(1.55) 



Cette equation peut se reecrire sous la forme 

p N+ i - 2p N + p N ^i c N+1 p N 



m 



(1.56) 
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Posons tj = t + je = t + j 



et definissons une fonction (p(t) par 



N + 1 

<p(tj) = €Pj, j = 0,...,N + l. 
II faut calculer lim (edet B) — lim (p(t). Or, ip(t) satisfait l'equation 

AT— »+oc N— >+oo 



cpitj + e) - 2ip{tj) + <p(tj - e) 



- -Cj+i 



me 



dV 
dt 2 



c(0 



m 



93 quand e — > 0, 



soit "i-nr + c(t)y? 
dt 2 



= 0. 



Les conditions initiales sont 

¥>(*o) 
d^(t ) 



ep -> 
Pi -Po 



df 



2 ci - 1 -> 1. 

m 



Le resultat final s'ecrit 



G(x,t;x ,t ) = 



m 



exp-[S, 



2irihip(t) h 



cl 



ou </?(i) est solution de l'equation differentielle 

d 2 y? 



m—- + V"(x c i(t))ip = avec < dy? 

i ^(to) 




1. 



(1.57) 



(1.58) 



C'est l'expression semi-classique. Elle correspond a un developpement asymptotique 
en h. 

Remarques : 



1. Cette approximation n'a de sens que si ip(t) 7^ 0! Si (p(t) = 0, il faut s'y 
prendre differemment. Le cas ip(t) = correspond a des systemes ou Taction 
n'a pas un minimum bien marque. 

2. Cette approximation ne peut etre definie que si la trajectoire classique existe. 
Ce n'est pas toujours le cas. En particulier les particules peuvent traverser 
des barrieres de potentiel en mecanique quantique. Dans ce cas, on peut 
encore faire un calcul approche du propagateur, mais il faut passer en temps 
imaginaire (voir chapitre suivant). 

3. Sous cette forme, la mecanique classique est a la mecanique quantique ce que 
l'optique geometrique est a l'optique ondulatoire. 



Chapitre 2 



Action euclidienne et temps imaginaire 



2.1 Definition 



Bien que revolution temporelle soit decrite par l'operateur e 
utile de connaitre le "propagateur euclidien" defini par 



l h H(t t )^ jj egt souven t 



G e (x,t;x ,t ) = (x e « 



x ), 



(2.1) 



ou r, r G IR. 

La terminologie "euclidien" vient du fait que Ton passe du propagateur standard au 
propagateur euclidien en remplagant it par r. Ce remplacement revient a passer de 
l'espace de Minkowski 



qui, au signe pres, est une norme euclidienne. Comme r = it, on parle aussi de 
"propagateur en temps imaginaire". 

Comme le propagateur est une fonction infiniment derivable pour t > to, le propa- 
gateur euclidien peut etre obtenu a partir du propagateur standard en remplagant 
t par —ir. En d'autres termes, on a : 



L'un des avantages du propagateur euclidien est que sa representation en termes 
d'integrale de chemin est mieux definie. Le propagateur euclidien est en effet donne 



X={t,X) =^ X =t —x 

a x = (r, x) =^ x = — t — x 



Ge{x,t;xq,t ) = G(x, -ir;x , -ir). 



(2.2) 
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par : 



G e (x,t;x ,t ) = 




r - t 
avec e = — — — . 
N 

Cette formule peut se demontrer directement en utilisant la formule de Trotter, ou 
en remplagant t par — ir dans l'expression de G(x,t; x ,t ) a l'aide d'une integrale 
de chemin. 

La difference fondamentale avec l'expression du propagateur standard a l'aide d'une 
integrale de chemin vient du fait que la quantite sous l'integrale est un nombre reel 
et positif. II peut done etre interprets comme une densite de probability sur les 
chemins. D'un point de vue mathematique, le statut de cet objet est mieux defini, 
et la convergence des integrates peut etre demontree moyennant certaines hypotheses 
sur V(x). 

Si Ton veut passer a la limite continue, la fonction qui intervient n'est plus Taction 
mais V action euclidienne definie par 

m— + V(x) 

ou L E (x) = m 9 ^- + V(x), le lagrangien euclidien, se confond avec l'hamiltonien dans 
le cas present. En termes de cette fonctionnelle, on a : 

G e (x,t;x ,t ) = [ Dxe-*^. 



Outre le fait que e'est un objet qui conduit a une integrale de chemin mieux definie, 
le propagateur euclidien apparait naturellement dans plusieurs contextes : 

• Mecanique statistique (fonction de partition) ; 

• Mouvement brownien ; 

• Effet tunnel. 
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2.2 Mecanique statistique quantique 



En mecanique statistique quantique, l'objet central est la fonction de partition 

Z = J2 e ~^ En = Tr exp(-/3if). 

n 

Considerons le cas d'une particule se deplagant dans un potentiel unidimensionnel. 
II est habituel de representer la trace dans une base d'etats propres \ip n ) : 

Z = Y t (<Pn\e- (}H \<Pn) 
n 

mais on peut aussi utiliser la base continue \x). En injectant la relation de fermeture 



on peut ecrire 



Z = / dx( if n \x ) ( x\e 

n J 



<Pn, 



X\e PH Wn) (<Pn\x] 



X 



Mais (x\e~ pH \x) = G E (x, (3%; x, 0) 



Z= dx Ge(x, @h; x, 0) 



ou encore 



Z = J dx G(x, — \(3h; x, 0) 



La fonction de partition a done une representation en termes d'integrale de chemin 

r x',pn 



= / d*< / 

J Jx',0 



Vxe~* SE[x] . 



(2.3) 



Exemple : Oscillateur harmonique. Puisque nous avons deja calcule le propagateur 
de Feynman de l'oscillateur harmonique, utilisons-le pour calculer sa fonction de 
partition. 



e /3hu, _ e -phuj i 



sm.{—\(3huS) 
cos(-iphu) = cosh(phu) 



2i 
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h Scl ~ 



avec xi = x f = x' . Mais cosh x — 1 = 2 sinh 

2mw 



x' 2 [cosh((3huj) - 1], 



,2 x 
2' 



h 



Le prefacteur s'ecrit par ailleurs 



hsmh((3huj) 

i 

mw \ 2 



sinh 2 f ^] x' 2 



2nh sin 



enne, on en deduit: 

dx'G(x', -iph; x', 0) 



h(/3hw)}) 



Par une integration gaussi- 



irh sinh(f3huj) 



2mu sinh 




C'est bien le resultat connu. En effet, 

+oo 



-f3hui/2 



n=0 



Z = 



e W2_ e -^/2 2 sinh (^) 
Limite classique : 

Dans l'expression de Z sous forme d'une integrate de chemin: 

r x',/3h 



dx' / Vxe-* SE W, 

Jx',0 



h apparait a deux endroits : dans la limite de l'integrale, et comme prefacteur de 
Paction. La limite h — ■> n'est done pas simplement une approximation du col. 
Quand H — > 0, les chemins qui s'eloignent de x' pour y revenir apres un temps 

d 2 x 

infiniment court sont penalises par le terme — -, et les chemins dominants restent 

dv 1 

proches de x'. Au premier ordre, on peut done remplacer V(x(t)) par V(x') 

rx' f3h 



x',0 



ou S%[x] est Taction euclidienne d'une particule libre. Mais le propagateur d'une 
particule libre a deja ete calcule (cf page 11) et est independant de x': 



G(x', —ifih; x\ 0) 



m 



m 



2n(3h z 



2.2. MECANIQUE STATISTIQUE QU ANTIQUE 



29 



Or. 



dpe 2m 



2nm 



Z r \ — 



2nh 



dpdq e" 



-PH(p,q) 



Au facteur 



1 1 



— pres, qui donne une mesure dans l'espace des phases, c'est la 



2nh h 
statistique classique de Boltzmann. 

Pour obtenir les corrections quantiques, il faut considerer, dans le propagateur 
G(x', t; x', 0), les chemins qui s'eloignent peu du point x'. Posons x(t) — x' + rj(r). 
Si t)(t) est petit, on peut faire un developpement de V : 



V{x' + r]) = V(x') + r]V'(x') + -r) 2 V"(x') + 0(r) 3 ) 



I V'jx'f 

= V{X) -2V^) + 2 V {X) 

-i /"* / V'(x') 2 
G(x',t;x\0) * exp-^ (v( X >) - 



V\x') 
V"(x>) 



+ 0{rf) 



At' 



V'(x') V'(x') 
V^^VHx 1 )' 



ou G osc est le propagateur de l'oscillateur harmonique dont la frequence est donnee 

11 
par -V"(x') = -muj 2 . 



x exp 



imw / V 



h sin ut V V 



(cos cut — 1) 



ou Ton a introduit les notations V = V(x'), V = V'(x'), V" = V"(x'). 

Remarque: V'(x') ^ car x — x' n'est pas la trajectoire classique. 

Comme on cherche les corrections a la limite U petit, on peut faire un developpement 
autour de t = (t — > —if3K). 
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sin cot 
cos out 
G(x',t;x',0) 



out — 



3+3 



6 



1 - 

exp 



ujH 2 ujH 4 

+ 



2 

-it 



24 

, V 

h V 2V" 



m 



2niht /-. _ u 2 t 2 

1 6 



x exp 



imuj 2 t 2 /2 ( V 



T777 -1 + 



uH 2 
12 



it V 
h 2V 



=> G(x',t;x',0) ■ 
On remplace t par —i/3h : 
=>• -i/3ft; x, 0) 



exp 



— iiV / m 



^(-i-^rl^) puisque ™ 2 =y" 

it 3 \/ /2 \ 



ft V 27ri/it 



1 + 



12m 



exp 



24hm 



exp(-W^/^(l-^-) (l + -2^- 



2tt/? 
I 1 -f 



12m 

I H7T „ / \ I 1 // 2 I '" f 2 lr 

\2m~ 



h\ 2np \ ' 24m 
Revenons a l'expression de la fonction de partition 



Z = J dx'G(x',-iph;x',0). 



Le propagateur G(x', — i(3h; x', 0) depend de x' via V = V(x'),V = V / (x'),V'" = 
^"(x'V Mais 



/+oo /"+oo r 

e-WJy"^)^' = _ / -/^'(x'K^VV) + e- pv ^V'{x') 
-oo •/— oo 



-oo 

f+OO 



! +oo 



/+OO 
-oo 

si le potentiel diverge a l'infini, e'est-a-dire si la particule est confinee dans une 
region de l'espace. Dans ce cas, on en deduit: 




C'est le debut du developpement de Wigner-Kirkwood. II contient les corrections 
d'ordre h 2 a la fonction de partition classique. 
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2.3 Mouvement brownien 



L'interpretation du propagateur comme somme sur tous les chemins possibles est 
intuitive mais purement formelle. Dans sa version euclidienne, le propagateur cor- 
respond effectivement a une somme sur les chemins possibles dans le cadre du mou- 
vement brownien. 

Pour decrire le mouvement d'une particule en suspension dans un liquide, Einstein 
a propose d'abandonner une description purement deterministe et a suggere une 
description probabiliste du mouvement : a chaque instant, la particule a une certaine 
probabilite d'effectuer un saut, et la description physique est basee sur la probabilite 
pour une particule d'aller d'un point a un autre en un temps donne. 

Dans le cas d'un mouvement unidimensionnel, l'hypothese la plus simple est d'admet- 
tre que la particule a une probabilite \ de faire un saut d'amplitude a pendant un 
temps r a droite ou a gauche. Si on definit la probabilite P(na, kr) pour la parti- 
cule de se trouver au point na a l'instant kr (n et k entiers), on a alors l'equation 
fondamentale : 

P(na, (k + l)r) = - [P((n - l)a, kr) + P((n + l)a, kr)} . (2.4) 
2 

Dans la limite ou a et r tendent vers 0, on peut faire apparaitre des derivees en 
retranchant de chaque cote P(na, kr) : 

P(na, (k + l)r) - P(na, kr) = ^[P((n - l)a, kr) - 2P{na, kr) + P((n + l)a, kr)} 

dP, . , a 2 d 2 P, , . 
r—(na,kr) = — — (na^kr). 

a 2 

Si a et r tendent vers de telle fagon que D = — tende vers une valeur finie, la 

2r 

limite continue de P satisfait l'equation differentielle : 

dP ^d 2 P 

W = D T^- <2 ' 5) 

C'est l'equation de la diffusion. Le coefficient D s'appelle coefficient de diffusion. 



NonnaUsaUon: f A *P( X ,t) = l. 
Condition initiale : lim P(x, t; xq, to) = 5(x — x ) 

t— >to 

(la particule est en xq a to)- 
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Cette equation est satisfaite par : 



P(x,t) 


1 

Vt 


x 2 
4Dt 








dP 

En effet, — 
dt 


1 . 


-3/2 e - 




-5/2^_ - 


4Di 


dP 

dx 


= r 1 ' 2 


—X 

e 

2Dt 


x 2 
4Dt 






d 2 P 

dx 2 


-l 

= 


-3/2 e - 


x 2 

ADt +t~ 


-5/2 £_ - 

4D 


x 2 


D 92P 

dx 2 


dP 
~dt' 











Avec la condition aux limites et la normalisation, on trouve 



On reconnait le propagateur euclidien dans lequel on a remplace m par Le 
propagateur euclidien correspond done a la probability de transition d'un point a 
un autre pour une particule classique brownienne. 

Le mouvement brownien permet de definir une mesure dans l'espace des chemins. 
Considerons en effet tous les chemins partant de x ,to et arrivant a x,t, et consi- 
derons N temps intermediaries tels que 

t <t ± < ■ ■ ■ <t N <t. 

La valeur moyenne d'une fonctionnelle de la forme 

F[x(ti), . . .,x(t N )] 

prise sur tous les chemins browniens allant de x , t a x, t est definie par 

(J 7 ) = J dx 1 J dx 2 ■ ■ ■ J dx N 

XP(x 1 t 1 , X t )P(x 2 t2, . . . P(xt, X N t N ) J r [x 1 , ...,X N ] 

On note cette valeur moyenne 

nXT 

dWF, 



J xn 



'So T 

ou dW definit la mesure de Wiener. Elle satisfait 

fXt 



pxt 

I dW = P(xt, x t ) (a verifier). 

J Xoto 
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Si on a une fonctionnelle plus compliquee du type 

tT 1 - ft drVlxM] 

la valeur moyenne est definie par une integrale de chemin 



v 



rxt ( N \ 

(xt,x to)= lim / dW exp | -e V[x(t k )] ) . 



Cette fonction P v (xt,x t ) est le propagateur euclidien d'un hamiltonien dont le 
potentiel est V(x). II satisfait done l'equation differentielle 

d d 2 

—P v (xt,x t ) = D— P v (xt,x t ) - V(x)P v (xt,x to) 
kmPy(xt, xoto) = 5(x — xo). 

Ce resultat est connu sous le nom de formule de Feynman-Kac. II se demontre 
exactement comme l'expression du propagateur en terme d'integrale de chemin : 
Pv(xt,x t ) = lim ... est precisement ce qu'on obtiendrait a partir du propaga- 

teur euclidien pour un potentiel V en utilisant la formule de Trotter. 
Remarques : 

2 

1. Le fait que ^- soit constant dans le passage a la limite continue montre que 

le long d'un chemin typique parcouru a vitesse constante, Ax oc V 'At. On 
Ax 

voit done que lim — — diverge =^ les chemins browniens ne sont pas derivables. 
At 

C'est l'origine des problemes rencontres dans la definition de l'integrale de 
chemin avec le potentiel vecteur + A(x i+1 )). 

2. Si V(x) est la probability que la particule soit detruite pendant une unite de 
temps, comme 

^11 (' - ^V-(x( ti ))) = ^exp (-<Ll*0 jy {x ( U ))\ , 

k=0 v 7 \ k=0 / 

le calcul precedent s'applique. La formule de Feynman-Kac a done ramene le 
probleme d'une integrale fonctionnelle a la solution d'une equation differen- 
tielle. 



2.4 L'effet tunnel 



Considerons le potentiel 

V(x) = 9 j(x*-a 2 ) 
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et essayons de calculer le propagateur 

G (a, +00; — a, — 00) = 



avec S[x] = / 

J — 



Vxe* s[x 



mx 2 — V(x) 



Si Ton essaie de faire une approximation semi-classique, il faut trouver une trafectoire 
satisfaisant l'equation 

mx = —V(x) 
x (—00) = —a 
x (+00) = a. 

Malheureusement, il est clair qu'il n'y a pas de trajectoire classique satisfaisant ces 
equations: pour etre en x = —a a t = —00, la particule doit etre au repos, et elle ne 
peut pas passer par dessus la barriere. Par consequent, l'approximation classique 
pour G (a, +00; —a, —00) n'est pas definie. 

II est neanmoins possible d'obtenir une approximation du meme type en considerant 
la version euclidienne : 



Ge (a, +00; —a, —00) 
S E [x] 



I 



Vxe~* SE[x] 



T 

2 



^mx 2 - V E (x) 



Le potentiel effectif pour ce mouvement est egal a Ve(x) = —V(x). 
V E (x) 




II doit done exister une trajectoire 
classique pour aller d'un maximum a 
l'autre. 



Cette trajectoire doit satisfaire 




-V'e(x) = V'(x) 
-a 



La fonction x(t) = atanh 



u(t - t c ) 



muo 



quelconque 



g 2 a 2 
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est solution du probleme. En effet, 



x = tanh 2 — — 



2 

x = xx = — ^x(x —a ) — — V (x). 

On peut done essayer de faire une approximation du col autour de cette solution 
pour obtenir une expression approchee du propagateur euclidien, et done du prop- 
agateur de Feynman. II se trouve qu'il y a encore une difficulte liee au fait qu'il 
y a une infinite de solutions puisque t c est arbitraire. Du coup, l'equivalent de la 
fonction (p(t) de l'approximation semi-classique est egale a 0. On ne peut done pas 
utiliser aveuglement la formule semi-classique. C'est le probleme des instantons (voir 
exercices). 



Chapitre 3 

Transitions de phase et integrale de chemin 



La theorie moderne des transitions de phase, et en particulier le calcul des exposants 
critiques, est basee sur l'utilisation de la theorie quantique des champs pour evaluer 
la fonction de partition. Le lien se fait en exprimant la fonction de partition comme 
une integrale de chemin pour un champ scalaire. Comme cette theorie des champs 
est la plus simple (par opposition aux theories de champs pour les fermions ou aux 
theories de jauge), il est desormais habituel de presenter la theorie des champs en 
commengant par cette application. 



3.1 Modele d'Ising: Introduction et generalites 

Les problemes rencontres en physique statistique sont souvent des modeles sur 
reseau, l'application la plus ancienne - mais toujours d'actualite! - etant le mag- 
netisme. Dans ce chapitre, nous allons voir comment on peut exprimer la fonction 
de partition d'un tel modele a l'aide d'une integrale de chemin sur le cas particulier 
du modele d'Ising. Dans ce modele, on represente le spin aux sites % d'un reseau par 
une variable <7; qui peut prendre les valeurs ±1, et l'energie d'une configuration est 
donnee par 



A? 
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L'aimantation sur site ( <7j ) = ^ l3E ^ (7i ^ 



Z 

• La fonction de correlation a deux points: ( <Ji<jj ) = — cr i i7je~^ E ^ CTl ^ . 

En general, il sufiit de resoudre ce probleme pour un champ magnetique uniforme 
(Bj = B independant de j) et dans la limite B — > 0, mais le calcul de Z pour B 
quelconque n'est pas plus difficile, et on en deduit aisement ( a { ) et ( a^aj ). En effet, 

dZ 



De meme, 



dB t 

W = 

d 2 Z 
OBidBj 

( CTiCTj ) 



^2 (3ate 



fo} 

1 <9Z ld(lnZ) 
/3ZaB~ ~ P dBi 



/3E(W) 



2 



1 1 dZ 

p^Z OBidBj 



II est done suffisant de connaitre Z({Bj}) pour connaitre toutes les fonctions de 
correlation. Z({Bj}) est appele fonction generatrice des fonctions de correlation. 

Enfin, si ( fTj ) ^ 0, la correlation entre les spins est mieux representee par la quan- 
tity : 

Gij = {(JiO-j) - (<Ji) (o-j) 

1 1 d 2 Z 1 1 dZ dZ 
SMt lj ~ P 2 Z dBidBj ~]PZ2dB i dBj 



Mais 



d 2 In Z d ( 1 dZ 



dBidBj dBj \ZdBi 



1 dZ dZ 1 d 2 Z 
+ 



Z 2 dBj dBi Z dBidBj 



1 d 2 In Z 
T,J 1 p 2 dB~dB~ 



s'appelle la fonction de correlation connexe, une terminologie qui vient du calcul 
perturbatif de cette quantite, et In Z s'appelle la fonction generatrice des fonc- 
tions de correlation connexes. Pour toutes ces quantites, Z({Bi}) est la quantite de 
reference. 



3.2. CHAMP MOYEN 
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Pour discuter la physique, il est plus intuitif de changer de variables et de decrire le 
systeme en termes des aimantations m ; = ( <7j ) plutot que des champs magnetiques 
B, L . Cela se fait a l'aide d'une transformation de Legendre qui fait passer de l'energie 

libre F({3; {Bi}) = — — In Z((3; {Bi}) au potentiel thermodynamique (ou de Gibbs) 



m {mi}) = F 09; {Bidmj})}) + £ B, 



rrii 



OF 

ou la relation Bi({rrij}) est donnee implicitement par rrii = —7775-- 
Les derivees de T sont egalement utiles : 



dr 

• — 

drrii 


dF 
^ dB~ 


dB j 
drrii 


„ x ^ dBj 


or 

drrii 


=~ m j 

= Bi 






d 2 r 


dBi 






drriidrrij 


drrij 







Mais si on considere comme fonction des rrij, eux-memes fonctions des B k , on a 
evidemment : 

dBi(mj(B k )) \ - dBi drrij 



^ dm,- 9i?fc 



e 2 r /?G, fc 



On en deduit que 

d 2 r 1 



drriidrrij f3 



3.2 Champ moyen 



La premiere chose a faire dans un modele de physique statistique est une approche 
de type champ moyen ou Ton suppose que les fluctuations des degres de liberte par 
rapport a leurs valeurs moyennes sont petites. Plus precisement, dans le cas du 
modele d'Ising, on definit 



5(Ji — (Ti — ((Ji) 
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et on neglige les termes d'ordre 2 en 5a : 



0i0 - = [5 Gi + ( <7j )) (5a j + (<Jj)) 

~ (a i )(a j ) + 5a i (a j ) + 5a j (a i ) 

= aiiaj) + a j (a i ) - (a i )(a j ) 

d'ou E = --^2j ij (a i {a j ) + a j {a i ) - (a^iaj)) - ^B^i 



3 



On s'est ramene a un systeme paramagnetique de spins decouples 



avec 



et 



Z,= ^ e /3B ^=2cosh(/5 J Bf) 



(Ti=±l 



. . 1 <91n Z 
mi = ( tTj ) = - = tanh 

13 dBi 



L'aimantation est done reliee au champ et a la temperature par l'equation de champ 
moyen 



rrii — tanh 
En champ nul, on a en particulier: 



P(B t + ^2 J ij m j) 



rrii = tanh | (3 Jij m j 



La configuration ou m; = pour tout % est toujours solution. Mais il doit exister 
une autre solution a basse temperature. Considerons pour etre concret le cas d'un 
reseau cubique en dimension D avec couplages J entre premiers voisins uniquement 
(J > 0). A T = 0, la configuration de plus basse energie correspond a wij = +1 
(ou rrii = —1)- Cherchons done a partir de quelle temperature il existe une solution 
uniforme 

rrii = m 

Elle doit satisfaire l'equation 



m = t&nh((3m2DJ). 



3.2. CHAMP MOYEN 
tanh() 

(32DJ- 
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m 

Cette equation a une solution non nulle si et seulement si (32DJ > 1 

=>• T < 2D J. 

II y a done une temperature T en dessous de laquelle le champ moyen predit 
l'existence d'une aimantation spontanee. 

Pour expliciter le contenu physique de la solution champ moyen, il est commode 
(mais pas indispensable) de calculer le potentiel thermodynamique. 



La relation champ moyen 



m, = tanh 



s inverse en 



Bi = — tanh 1 m 8 — J, 

" 3 

L'energie libre est donnee par : 

F(f3;{B,}) = -ijVln2-i^lncosh 

4l> 



^,3 



T(P;{m t }) = F{p-{B l {m 3 )}) + Y,B l m l 

i 

= --iVln 2 - - In cosh tanh" 1 (m^) + - ^ J > 



P 



1,3 



+ — mj tanh 1 rrii — m; Jy-m.,- . 



Or, cosh (tanh _1 x) = (1 - x 2 y 1/2 . En effet, 

y = tanh" 1 x ^ 1 — x 2 = 1 — tanh 2 y = 



cosh 2 y — sinh 2 y 
cosh 2 y 



cosh 2 y 
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1 1 + x 

De plus, tanh _1 (:r) = - In . En effet, si y = tanh -1 x, 

F ' w 2 1 -x ' y 

1 1 + tanh y 1 cosh y + sinh y 1 2e y 1 , 2 , 

- In — = - In — — = - In = - In e y ) = y 

2 1 - tanh y 2 cosh y - smh y 2 2e~^ 2 v ' 



T{f3-{mi}) = ~Y^7n i J ij m j -^hi2 

1 



Susceptibilite uniforme : 

Si on s'interesse aux solutions uniformes, to = et pour un reseau cubique en 
dimension D avec J > entre premiers voisins, il vient : 

r(/3;m) = -TV JDm 2 - ^ In 2 







AT 



H — - {(1 + to) ln(l + to) + (1 - to) ln(l - to)} 
2p 

Le potentiel de Gibbs par site 7 = T/N est donne, pour to petit, par : 



m 2 T 
7 (/3;to) = -Tln2 + — (T-2JL>) + — m 4 

I 

, /-, x m2 m3 ^/ 4\ 

puisque m(l + TOj = to — + — — hC(m J 



d'ou 



TO 2 T 

7 (/3; to) = — T In 2 + — (T — T ) + -m 4 . 



7 / 


>■ 




T>T 




TO 



7 ' 


x7 , 


T<T \ 


V/ 




TO 



II y a done une transition de phase a T c puisque le minimum se deplace de to = a 
une valeur finie (deux fois degeneree). 
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Par definition, la susceptibilite magnetique est la derivee de l'aimantation par rap- 
port au champ en champ mil. Son inverse est done donne par: 



X 



dB 
dm 



<9 2 7 



m=0 



dm 2 



T-T 



m=0 



soit 



X 



(T-Tq)- 1 . 



De fagon generale, des arguments de lois d'echelle predisent que la susceptibilite doit 
se comporter comme x oc (T — T c )~ 7 au voisinage de la temperature critique, ou 7 
est un exposant critique. A l'approximation du champ moyen, l'exposant critique 
7 = 1- 



Fonction de correlation : 



Comme l'approximation champ moyen consiste a negliger Sa^Sa j, elle ne peut etre 
precise que dans un domaine de temperature ou ( ba^aj ) est petit (condition neces- 
saire). Calculons done 

( 5a i 5a j ) = ( a,^ ) - ( a, t ) ( a j ) = G tj . 
Ce calcul se fait aisement a l'aide de T puisque 



(3 dmidrrij ' 



soit 



En effet, 



i(l + x)ln(l + x)] = ln(l + x) + l 

d 2 1 
-[(l + x)\n(l + x)} = 



d 2 x 1 + x 

Mais dans la phase haute temperature, rrii = =>- T (G" 1 )^. = — + T5ij. 
Toujours dans le cas d'un reseau cubique en dimension D avec couplages entre 
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premiers voisins, on peut introduire la transformee de Fourier 

D 



J(q) = Yl eim ~ rl)j v = 2J Y, C0S( i° 

i 

T{G~ 1 ){q) = J^e if - (fi -^ ) T(G- 1 ) t 

i 

D 



a=l 
■ij 



a=l 



G(q) = 



T 



T - 2J Ef=l C0S Qa 



d'ou 



dq- 



(2tt) d J zb T-2JEf =1 cosg a 

Cette integrate est dominee par les petites valeurs de q (le denominateur est plus 
grand) 

* Gl3 ~ (2^ J zb ^T-T + Jq^ 
Considerons le cas D=3, et posons f = fj — fj. 



T 1 f +co q 2 dq 



J(27T) 2 i £- 2 + <7 2 io 



J 

sin^e- igrcos9 



avec £" 2 = 



T I f +0 ° a 2 da f 1 



T 1 f +00 q 2 dq 2 sin qr 



Til A 4 " 00 gdgsin qr 



J(27r) 2 r./-oc ^ 2 + g 2 
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., ,. ; +0 ° qe iqr dq 

Mais " ' 



/ gdgsin qr f 



oo 



r i e -/€ 



= 7T — ■ = 7TC 

Ce dernier resultat s'obtient en passant dans le plan complexe, et en fermant le 
contour par le haut pour englober le pole a q = i^ 1 . 

Ainsi, G(r) = e~ r/f 

£ = \J^f~~^f~ l° n g ueur de correlation. 

Les lois d'echelle predisent que la longueur de correlation diverge comme (T — T c )~ u , 
ou v est un nouvel exposant critique. On a demontre que l'exposant critique v — \ 
a l'approximation du champ moyen. 

Loin de T c , ^ est petit et Gij est tres petit. Dans ce cas, le champ moyen est 
une bonne approximation. Mais pres de T c , les spins sont correles sur de grands 
distances, G^ n'est pas petit, et le champ moyen ne peut pas etre bon a priori. 



3.3 Modele d'Ising et integrale de chemin 



Le champ moyen tel qu'il a ete formule dans la section precedente ne donne pas 
naturellement naissance a une methode permettant d'inclure des corrections. Cela 
peut etre fait de fagon systematique en reformulant le probleme comme une integrale 
de chemin. 

Essayons d'exprimer Z = ^ e~ fSE a l'aide d'une integrale multiple. On peut ecrire 

o-pE _ ±*<t(/3J)*+I3*B<t 



ou on a introduit les notations 

a 



J 



(<7i, . . . ,a N ) 



( Jn J 



11 -Jl2 
J21 J22 



J 



k B = (B 1 ,...,B 



N 
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II est commode d'inclure la temperature inverse f3 dans les matrices en definissant : 

A = pj 
H = (3B 

z = Y^^ taAa+tHcj ■ 

L'idee de base est d'utiliser l'integrale gaussienne pour ramener cette expression a 
une expression lineaire en a, ce qui permettra d'effectuer simplement la somme sur 
Or, d'apres la formule de base de l'integration gaussienne, 



i r N 



1=1 



avec 



Z(0) = (2^)^-^= = (2n)»'*Vd*A 

=>• Z = ■ — . / dxi . . . dx N e~2 

(2n) N / 2 Vd^A J 

(on a utilise t aH =* Ha). 

Effectuons maintenant le changement de variable (p = x + H =>- x = <p — H : 

= ~\ WV + \ 'HA-\ + \ *<pA- x H - \ *HA- X H 
= ~ t ipA- 1 i f + t HA- 1 p-- t HA- 1 H 

La somme sur {a.i\ peut maintenant etre effectuee tres simplement 



<T 1 =±1 

cro=±l 



i <Ti=±l 



= ( J (<••• • <• '-■') 

i 

= 2 JV J]cosh(^). 
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On ecrit cette expression exp(— U(ip)) avec U(ip) = — lncosh(^j) — A^n 2. En 
remplagant A et H par leurs valeurs en fonction de J et B, on trouve finalement 

Cette expression est semblable a l'expression du propagateur euclidien, ipi jouant le 
role de x^. La difference majeure vient du fait que i prend des valeurs sur un reseau 
de dimension D, alors que dans le cas du propagateur, % prend des valeurs de 1 a A, 
done sur un reseau de dimension 1. Autrement dit, la forme generale de la fonction 
de partition a la limite continue est 

Z oc J Vcpe- SM 

ou Taction est une integrale multiple : 

S[<p] = [ d d fL[ip}. 



Champ moyen 

L'approximation la plus simple consiste a faire une approximation du col sur cette 
integrale multiple. L'exposant s'ecrit 

hi hi i 

OS 1 x 

Notons (pi la solution de ce systeme. L'approximation d'ordre pour Z s'ecrit done 

e -*M (3.1) 



(2tt) n / 2 Vdet [3 J 
A une constante pres, on a done : 

\nZ = ~ 'BJ^B - l (pJ- V +' BJ~ x {p + U((p) 

soit 

In Z = -| Y.B,J,/B, - YQ Y^rrl./rj + £ B * J ^i + U(if) 

hi hi hi 
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Or, l'aimantation est donnee par 



rm = (<Ti) 



D'apres l'expression de In Z, il vient: 



1 <9(ln Z) 



dBi 2 ~ J 2 

3 3 



Sous forme vectorielle, on a done: 



ou encore 



Mais au point col, on a: 



Ainsi, 



et 



B = —Jm + —Lp 



— J 1 ip — J 1 5 = tanhy9, 



m = tanh tp 



B = — Jm + — tanh 1 m 



Bi = — tanh 1 — Jij m j 



On retrouve bien les equations de champ moyen. 

Pour aller plus loin, la premiere etape consiste a calculer le determinant de la matrice 
dont les coefficients definissent la forme quadratique correspondant aux corrections 
du second ordre dans un developpement de S autour de (p. Ceci peut se faire 
directement sur la forme discrete de S, mais il est instructif de prendre la limite 
continue pour plusieurs raisons : 



• D'apres l'analyse de la solution champ moyen, ce sont les fluctuations a q petit 
qui dominent pres de T c . Une approximation continue est done justifiee. 
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• Le modele continu a des proprietes universelles : on peut demontrer que les 
exposants critiques ne dependent que de considerations generates (symetries, 
nombre de composantes du parametre d'ordre). En etudiant directement le 
modele continu, on etudie directement toute une classe de modeles. 



• C'est sous la forme continue que le lien avec la theorie quantique des champs 
deviendra evident. 



II faut done determiner la version continue, en champ nul, de S donnee par 



La constante est sans importance, et on s'interesse a une region ou ipi est petit (on 
considere les petites fluctuations autour de </?, et <p = pour T > T ). 




i 



In cosh x 




2 12 



In cosh ipi 



2~XX + 12 XX 



Par ailleurs, 



9192 i,j 



i{qi-q)n \^ „i(9-92)r 3 

N ^ 




3 



E ( J ~ _1 ) (?) 
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Mais 



D 



J(q) = 2 J ^ cos q a 



a=l 
D 



1 



T - Jg 2 
1 



T - Jg 2 T \ 1 _ 

\ To 

1 / Jq 2 
~ — 1 + — 
To T„ 



»J ° 9 ° 9 

J "^T* -^0 

1,3 hJ 



T 
1 T 
27^ 



soit 



1 JT 

Si on definit fa = \ I -^2 W 



1 



T 2 

— ^- </>j, cela s'ecrit 
JT r ' 



2 JT 



(T-T )J2 



1 T 4 



12 J 2 T 2 



5>f- 



On peut remplacer la somme par une integrate 





"1 




^ dr 


2 


V0 



2 + ^ O (T)0 2 + ^ O (T)0 4 



avec 



r (T) = r (T-T ), f 



2^ 



JT 



A un facteur multiplicatif sans importance pres, la fonction de partition s'ecrit: 



Z = / T>0e 



-S[4>] 
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C'est le modele dit de Landau- Ginsburg 1 . 



3.4 Modele de Landau-Ginsburg 



Le modele de Landau-Ginsburg est defini par la fonction de partition 

Z = J V(j>e~ m 



avec 



sm = J 



dr 



V0 



2 r 



no 4 
4 ^ 



- B{r)<t>{r) 



(3.2) 



ou r = r (T — T ), et ou on a ajoute un terme de source J dfB(r)(j)(r) a Taction 
par commodite 2 . 

Lorsque <fi est un scalaire, comme dans le cas du modele d'Ising, ce modele decrit 
de fagon generale les transitions de phase du second ordre dans les systemes ou le 
parametre d'ordre a une seule composante. Ce modele peut etre generalise a des 
parametres d'ordre plus compliques avec des consequences non triviales. 

Si uo = 0, on parle de modele gaussien. Ce modele n'est defini que pour r > 
(i.e. T>T ). 

Pour se familiariser avec la formulation continue, il est utile de refaire l'equivalent 
du champ moyen. Dans ce contexte, cette approximation s'appelle V approximation 
de Landau. 

Cette approximation consiste a chercher le champ qui minimise Sty], et qui doit 
done satisfaire la condition: 



avec 





"1 




j dr 


2 


V0 



50(f) 



1 







2 r o 2 + ^0 4 -5(f)0(r) 



1 Certains auteurs ecrivent le terme en d> 4 avec un coefficient — 7 

4! 

2 Le terme de source n'est pas vraiment un champ magnetique. C'est un intermediaire de calcul 
qui permet de calculer les fonctions de correlation de la theorie. De m§me, Bi n'est pas un champ 
physique. II permet simplement de calculer (<7j<7j ). 
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L'equation s'obtient par une generalisation triviale du cas unidimensionnel : 

SF[f] 



F[f] = JdxG(f(x))-. 
F[f] = J dxG(f(x)) 



Sf(y) 
SF[f] 
Sf(y) 



J dxG' (f(x))5(x-y) = G> (f(y)) 

dxG' (f(x)) 5>(x -y) = —G" (f(y)) f"(y) 



Ainsi 



soit 



50(f) 



= -A0 + r o + m o 3 - B(r) = 



-A0 + r o + w o 3 = B(r) 



Soit 0o la solution de cette equation. L'approximation de Landau consiste a poser 
Z = exp(— S'[$o])- Essayons de determiner le potentiel de Gibbs a cette approxima- 
tion. Dans ce contexte, l'aimantation est definie par 

^. . 5 In Z\B\ 5S[B] ^ 
M(r) = 1 J = - 1 J = $ (f) 



5B(f) 



SB(r) 



(le facteur -| a ete inclus dans le champ lorsqu'on a introduit le terme de source 
c'est la convention habituelle 3 ) 

T[M(f)] = F [B[M(r)]\ + J dfB(r)M(f). 

Mais a cette approximation F = S[& ] 



2 + I ro M 2 + — M 4 
2° 4 



Cette equation est la version continue de l'equation donnant T(f3, {m,j}) dans le cas 
du modele sur reseau. 

Dans le cas uniforme, on peut ecrire le potentiel de Gibbs par unite de volume sous 
la forme: 

7 = — = -r m 2 + -^-m 4 (champ moyen). 





"1 




j df 




VM 




2 





3 Cette relation est a rapprocher de rrii = taring. Comme tanh0i = 4>i — on peut montrer 
que cette difference est sans importance pour les quantites que nous calculons. En realite, c'est 
clair puisque r[M(r)] est bien la version continue de {m^}). 
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Pour aller a l'ordre suivant, il faut faire un developpement de S[<f>] au second ordre 
autour de O ; 



S[(f) + rj] = S[<f> ] + J drr)(r) 



6S[<j>] 



5(f){f) 



-<t>0 



S 2 S[4>} 



50(fi)(50(f 2 ) 



et il faut calculer 

J Vrj exp 

Ecrivons 



1 1 dPi I d ^ (r ' Mr ' 2) mmm 



<t>=<t>0 





Si[<t>]+S 2 [<j>] + S 3 [<j>]- 1 


SM = 






2 




v<t> 




5 2 [0] = 




4? 




= 


/*7 


<t>\ 





Le dernier terme ne contribue pas a la derivee seconde. Calculons les derivees 
secondes des autres termes. 

Pour simplifier, supposons qu'on soit a une dimension. 



6SM 
S<f>(xi) 

50(xi)50(x 2 ) 



dx 4>'(x)6'(x — xi) 

dx 5'(x — xi)5'{x — x 2 ) 



dxi / dx 2 r)(xi)r)(x 2 ) 



5<j)(x 1 )S(j)(x2) 



dxr]'(x) 2 
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5S 2 



5<p(xi) 
S 2 S 2 [<j)} 



5(f)(x 1 )5(f)(x2) 
J ctoi J dx 2 r)(x 1 )r)(x 2 ) 



/o / dx (f)(x)5(x — xi) 
r J dx S(x — xi)S(x — x 2 ) 
5 2 S 2 [<P\ 



5<f)(xi)8(j)(x 2 ) 



-<t>0 



J dxi J dx 2 rj(xi)rj(x 2 ) 



5^x l )5^x 2 )=^S dXri{xf - 
^jdx^x-x^x-x.) 

= 3u Jdx^Mxf. 



s 2 sm 



54>(xi)5(f>(x2) 



Finalement, il faut calculer 

Vrj exp 



\ J d r(-r]Ari + r ri 2 + 3u (j)lri 2 ) 



Supposons toujours O uniforme. C'est l'analogue en trois dimensions du probleme 
de l'oscillateur harmonique avec des conditions aux bords non limitees. 

Pour calculer cette integrale fonctionnelle, il faut en principe considerer une version 
discrete, calculer l'integrale multiple, puis prendre la limite continue. II serait plus 
logique de faire ce calcul avant de prendre la limite continue du modele sur reseau, 
ce qui est bien sur faisable. 

Neanmoins, la forme precise du prefacteur n'a aucune signification physique (on ne 
change pas la physique en multipliant la fonction de partition par une constante), et 
il nous suffit de connaitre la dependance fonctionnelle du resultat sur les parametres 
r , u et 0o- 

Or, si on passe sur le reseau, elle est donnee par (det A) -1 / 2 , ou A est la matrice 
decrivant l'exposant. Mais det A = le produit des valeurs propres. II suffit done de 
connaitre les valeurs propres de l'operateur 

A = -A + r + 3w $o- 

Or, cet operateur est diagonalise par transformee de fourier, la valeur propre associee 
a q etant q 2 + r + 3m o 0o- En effet, 



Ve^'" = iqe 1 ^ 
Ae if " = igVe if " 

det A = Y[ (<? 2 + + 3m o 0o) • 



-q 2 e iq ^ 
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On en deduit 

In Z = -S(0 O ) - l - In (q 2 + r + 3u $) . 
<? 

On a ainsi calcule la premiere correction a l'energie libre 4 . 

Pour faire le lien avec l'approximation de Landau, il est utile de calculer le potentiel 
thermodynamique au meme ordre. Commengons par exprimer l'energie libre par 
site, toujours avec 0o uniforme: 



df 

m = — 



OB ' dB"Z ^ 

=o 



d<j) 1 f dq 1 

x 6u (p - c 



OB ^2J (2n) D q 2 + r + 3« $o 
= 0o — $ m avec $ m = 4>q — m. 
0o = m + 5m. 



Le potentiel de Gibbs par site 7 est donne par: 



7 = / + Bra = -r (m + 5m) 2 + -u (m + 5m) 4 



1 f 

-B5m + - J j^y5 ln (l 2 + r o + 3u (m + 5m) 2 ) 

= -r m 2 H — ^m 4 + 5m(r m + ur,m 3 — B) 
9. A .. ' 



-0(O(5m)) 



soit 



7 



-r m + — m +- J D ln(g + r + 3w ™ ) 

" ^ ' " ^ ' 

Approximation Premier terme du developpe- 

de Landau ment en nombre de boucles 



4 Ce choix correspond a un choix de mesure J T)r\ tel que le prefacteur soit egal a 1. 
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-i <9 2 7 

Pour etudier le systeme a cet ordre, il suffit de calculer x = ^ — o 

dm 2 

dq 1 



m=0 



— — = r m + u m 3 + 3-u m 
dm 



(2tt) d q 2 + r + 3w m2 



<9 2 7 



(9m 2 



r + 3u 



m=0 



dq 



(27r) D q 2 + r (T) 



La temperature critique T c est donnee par x 1 — 0, soit 

dq 1 



f (T c -T ) + 3 Mo J 



soit 



?~0 J 



)^g 2 + fo(T c -T ) 
3m f dq 1 



= 0, 



(27r)^V + fo(T c -To)- 



Cette approximation n'a de justification que si T c — T est petit, auquel cas on peut 
negliger f(T c — T ) au denominateur, ce qui donne 



T C = T 



3tin 



f 7 (27r) B g 



d<f 1 

D 72 



La temperature critique est abaissee par les fluctuations. 

Etudions par ailleurs le comportement de x~ X P r es de T c . 

dq 1 



X 



r (T - T ) + 3u 



(2n) D q 2 + f (T-T ) 
f (T-T c )+f {Tc-T^ +3u °Jj^ 



- 3uo f-^2_ i 

U J (27r) D q 2 




f (T-T c ) + 3 Mo 



d<f 



(2tt) d V^ 2 + X" 1 9 2 



Mais 



d"ou 



q 2 + X- 1 q 2 q 2 (l + ^) 



q 



in 



x 



X 



-1 



X-^foCT-Tj-S^oX- 1 



q , f 
dq 1 



7 (2n) D q 4 ' 

Le comportement de l'integrale depend de la dimensionnalite. 
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D > 4 : l'integrale converge. Posons 



C = 3u 



dq 1 



(2tt) d g 4 

L'exposant critique 7 est toujours egal a 1. On peut demontrer que cela 
reste vrai si on pousse le developpement aux ordres suivants. On resume en 
general cette situation en disant que le champ moyen est exact en dimension 
superieure a 4. En fait, T c est modifiee par les fluctuations, mais pas les 
exposants critiques. 

D < 4 : cette integrale diverge. Si on veut etudier le comportement tout pres de T c , 
on ne peut pas se limiter a cet ordre. La technique appropriee, le groupe de 
renormalisation, sort du cadre de ce cours. 



Si Ton se limite a des temperatures pas trop proches de T c , cette approximation 
reste neanmoins valable. 

Pour determiner la limite de validite, revenons a l'equation 

*" = fo(r - Tc) + 3l, °/(lF(?TFT-?) 

= f ° (T - T ' ) + 3 "°/(^(?(?T^)) 

Dans cette integrale, il faut traiter separement q < a/x _1 et q > a/x _1 

r\Jx q D ~ l dq 

Jo ( 27 T 



q> ^'- h K //-r(2W 



L'approximation reste valable tant que le terme correctif est beaucoup plus petit 
que le premier, soit tant que: 

u (f (T-T c ))^ 2 <l. 
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C'est le critere de Ginsburg. 

Remarque: l'analogie avec la theorie quantique des champs est evidente. En theorie 
des champs, c'est directement le champ qui a une signification. On n'a pas besoin 
de faire de transformation de Legendre (masse m l/£ longueur de correlation). 



Chapitre 4 



Integrate de chemin et probleme a A^-corps 



4.1 Introduction 



L'objet central du probleme a iV-corps est de calculer la fonction de partition 



n\ e 



-f5H 



n 



Pour obtenir une representation en termes d'integrale de chemin, il sufiit, comme 
dans le cas a une particule, d'utiliser la formule de Trotter et d'introduire une base 
continue a chaque etape. Le choix qui peut sembler le plus naturel consiste a utiliser 
les fonctions symetrisees ou antisymetrisees (bosons ou fermions) construites a partir 
des fonctions propres de Xi et pi. Cette repesentation s'avere utile dans certains cas si 
l'on veut essayer d'evaluer Z numeriquement en utilisant par exemple une methode 
de type Monte Carlo. On peut montrer que Ton arrive a une expression du type 

Z = j^/^£ P f dx 1 ...dx N J Xl( i3) = XPl(0 ) = X1 Vx 1 ...Vx N 



AO) = x N 



x exp 



pour un hamiltonien du type 



dr ^™^y + i Ev(I((T) _ Ij(r)) 

i=l ^ ' iy^j 



H 



N 

E 



Pt_ 
2m 



1 i^j 



La somme sur P est la somme sur toutes les permutations des indices % — 1, . . . , iV, 
£ = +1 pour les bosons et —1 pour les fermions, et P est la parite de la permutation. 
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Cette representation ne se prete pas simplement a des developpements analytiques 
a cause de la somme sur les permutations. Cette difficulty vient directement de 
l'indiscernabilite des particules qui oblige a travailler avec des fonctions totalement 
symetriques ou totalement antisymetriques. 

Pour arriver a une expression donnant lieu a des developpements analytiques in- 
teressants, il est naturel d'essayer d'utiliser la base de Fock. Mais la base de Fock 
est definie a l'aide de variables discretes : les nombres d'occupation. Pour obtenir 
une representation en termes d'integrale de chemin, il est utile d'introduire une base 
continue dans l'espace de Fock, la base des etats coherents. 



4.2 Etats coherents : bosons 



La construction des etats coherents ne se fait pas de la meme fagon pour les bosons 
et les fermions. Dans cette section, on s'interesse aux bosons. Commengons par le 
cas le plus simple, celui d'un oscillateur harmonique, et rappelons les proprietes des 
etats coherents. 

L'oscillateur harmonique est defini par l'hamiltonien 

p 2 muj 2 9 
H = — H q . 

II est possible de reecrire cet hamiltonien a l'aide d'operateurs a, a + definis par : 

1 



a = 



1 


Imuj 




[V n 


i 


Imuj 


7! 


LV n 



q + i 



q - l 



1 



On verifie que [a, a + ] = 1. 
Demonstration : 



[a, 



A l'aide de ces operateurs, le hamiltonien s'ecrit : 

H = hu ( a + a + - 
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Demonstration 



p = 1 Jme*(a + -a) 



p 2 m 2 u) 2 2 muifi , 2 +2 + _■_ N 

1 q = [a +a —aa —ara) 

2m 2 H Am V ; 

muj 2 k . o ,2 , , s 

H (a +a + aa + a a) 

2 2mu) 

----- iKiS [ (I ' (I ■ - J . 



• Les valeurs propres de a + a sont positives ou nulles car si a + a\(p n ) = a n \(p n ), 

(ip n \ a + a \ip n ) =a n = ||a|v? n )|| 2 > 0. 

• Si \ip n ) est vecteur propre avec a n comme valeur propre, a\ip n ) est vecteur 
propre avec a n — 1 comme valeur propre. 

Demonstration : 



a + aa\ip n ) = (— 1 + aa + )a\ip n ) 
= -a\ip n ) + aa n \ip r , 



Pour ne pas generer d'etats propres de valeur propre negative, il faut que le spectre 
soit entier pour que l'application successive de a annule l'etat de depart. On appelle 
|0 ) l'etat de valeur propre pour l'operateur a + a. 

Definition: Les etats coherents sont les etats propres de l'operateur a. 

Demontrons que e za+ |0) est un etat propre de a, et calculons sa valeur propre. 



ae za 



00 r 
z 



?1=0 



Or, 

[a, (a + ) n ) = [a, (a + )"- 1 ] a + + (a + )™- 1 = n(a+)"- 1 
(se demontre par recurrence). 

a(a + ) n \0) = (a + ) n a|0)+n(a + ) n ~ 1 |0) 
=0 

ae 2a+ |0) = V- r T (a + ) n - 1 |0) = ze za+ W 

n=l v ' 



62 CHAPITRE 4. INTEGRALE DE CHEMIN ET PROBLEME A N-CORPS 



(dans la somme, n = donne car a|0) = 0). Ainsi, \z\ 
propre de a avec la valeur propre z. Par ailleurs, 



= e za 10 ) est vecteur 



a\z } = z\z 



(z\a 



z\z . 



Relation de fermeture: A l'aide de ces etats coherents, qui constituent une base 
surcomplete, on peut demontrer la relation de fermeture suivante : 



dz* dz 



ou la "mesure' : 



2ivr 

dz*dz . f d(Re z) d(Im z) 
est dennie par / . 

2l7T J 7T 



* z \z)(z\ = 1 



Demonstration : soit \n) la base orthonormee des etats propres de a + a ayant la 
valeur propre n. On a 



|n) 



10). 



En effet, 



(0|a n (a + ) n |0) = 



(0|a™- 1 (a(a + ) n )|0) 
(0|a"- 1 (n(a + ) ri - 1 + (a+)"a) |0 

n(0|a n - 1 (a + ) n - 1 |0) 



n! par recurrence. 



On peut done ecrire 



et 



n=0 



+°° 



Ainsi 



/ 



dz* dz 
2m 



e- z * z \z)(z\ 



m=0 



d(Rez)d(Im z) ^" , w 

— — -e } j —j=—j=\n){m\. 



nl Vml 



Passons en coordonnees polaires : z = re 



i0 



dz*dz ..„,,, /"rdrdfl r « e i0v m e- i9m 



2i7T 



\z){z\ 



IT 



n\ym\ 



n){m\. 



Mais 



1 d^e i ^- m ) = 2 7 r ( 5 m , n , 
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d'ou 

dz*dz „.„, , , , /" , v-e- r V 2n 



/d.Z* dz * /" g-r- a r 2n 
e~ z * z \z) (z\ = / 2rdr > ; — 
2i7r 1 M 1 J ^ n\ 1 A 1 



II reste a calculer 

/„ = / 2/--" +i <> '-(I/- 



dw w n e u (u = r 



2^ 



/■+oo /■ 

/ duu n e~ u = - du (-e~ u ) nu* 1 - 1 + [-e- u u n ]+°° 

Jo J v v ' 

=0 



/ 



= nl n -i = n\ 
dz* dz 



2m 



e z * z \z) (z\ = \n)(n\ = 1 



c.q.f.d. 



Produit scalaire: Les etats coherents ne sont pas orthogonaux. Leur produit 
scalaire est donne par (z 2 z 1 ) = e z i Zl . 



Demonstration 



— ym! v — v — ' v Til 



(4) m / , \ zi 

8„ 

E 



n 



Fonction de partition : 

Considerons desormais un systeme decrit par un hamiltonien qui est un polynome 
en a + , a. Le but est de calculer sa fonction de partition exprimee sous la forme: 

Z = Tr(e-? H ) = J ^^e- zz *(z\e^ H \z). 

On est done amene a calculer la valeur moyenne d'operateurs entre des etats co- 
herents. Or, ce calcul n'est simple que si les operateurs sont dans l'ordre normal, 
e'est-a-dire si, dans chaque terme du polynome, les operateurs a sont a la droite des 
operateurs a + . Dans ce cas, on a en effet: 



z 2 \(a + ) m a n \ Zl } = (z* 2 ) m z™ (z 2 \ Zl ). 
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Par contre, si ce n'est pas le cas, il faut remettre les termes dans l'ordre normal 
en utilisant les relations de commutation. Mais meme si le hamiltonien est dans 
l'ordre normal, son exponentielle ne Test pas puisqu'elle fait intervenir les produits 
des differents termes du polynome. 

Pour resoudre ce probleme, on utilise une variante de la formule de Trotter: 

e"^= lim (t-^ 



N^+oo V A" 

et on introduit une relation de fermeture entre les differents facteurs, ce qui conduit 
a l'expression: 

Z = lim [*^ e -* /^-i^-ie-^ 
w-»+oo J 2ixi J 2m 

^-«<,|i-^|, w _ 1 >...<* 1 |i-f|*> (4.i) 

On est done ramene a revaluation de valeurs moyennes du type 

/ Ml P R \ \ 
{z 2 \l-—\ Zl ) 

Or, si H est dans l'ordre normal, on a: 

(z 2 \H(a + ,a)\z 1 ) = H(z*,z 1 )(z 2 \z 1 ). 



I 



Ainsi, 



lim f dzdz * f dz i dz *i f 
n^+oo J 2ni J 2ni "J 



avec 



2vri 

xe -| [H{z* ,z N - 1 )+-+H(z%,z 1 )+H(z* 1 ,z)\ 

lim ff f dz k dz t e -s(zt, Zl ,...,z* N ,z N ) 
Af^+oo-Ll / 27ri 

k=l J 



N ( P * \ 

fc=i ^ ' 



et z = Zjv- Mais 



*/ \ P * f Zk ~ z k-l\ 



(3 d 

— z*(t)—z(t) quand A" — > +oo. 

AT v ' dr v ' 
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Finalement, on peut introduire la notation suivante : 



Z= I Vz*Vz exp[-S[z*,z]] 

lz(J3)=z(0) 



avec 



S[z*(t),z(t)}= f 
Jo 



z*(r)^-z(r) + H(z*(r),z(r)) 



dr. 



On a done exprime la fonction de partition comme une integrate fonctionnelle. La 
mesure est definie par 

Vz* Vz = f[ dZk dZ * k = f[ ^— d( - Im 



2m 

k=l k=l 



7T 



La notation n^Li ^ifi vient d'une autre fagon de voir les choses. On peut definir 
une integrate de chemin sur l'espace des phases dq dp, puis faire un changement de 
variables 

z=^=(Q + ip), z* = -^=(q-ip) 
d(Re z) d(Im z) dpdq 

7T 2li 

Or, si on considere z(p,q) et z*(p,q) comme un changement de variables, on a : 



dz dz* = dp dq x 



= dpdq x 
dp dg dz dz* 



dz dz 

dp dq 

dz* dz* 

dp dq 

i 1 
-i 1 



dpdg x i 



2tt 



2m 



„ ni T ,.r>-- t . f d(Re z) d(Im z) 
Cette manipulation est purement formelle ! La definition est bien / . 

J 7T 

De meme que l'integrale gaussienne est l'element de base pour le calcul des integrates 
de chemin de Feynman, la formule de base pour l'integrale de chemin des etats 
coherents est : 

i=i 



ou H est une matrice hermitienne positive (voir appendice mathematique). 
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Exemple : Calculons la fonction de partition de l'oscillateur harmonique. 



Z = lim 

7V->+oo 



JV 

n/ 

k=i J 



dzt dz, 



fc uz V zt A a*i ( e -N^) N 
2m v ; 



avec 



A = 




et a = 1 — —hu 
N 

Un developpement du determinant suivant la premiere ligne conduit a: 

det A=l + (-l) JV - 1 (-a)(-a) JV " 1 = 1 - a N . 



Mais 



lim ( 1 + 



N 



Ainsi, lim det A = 1 - e~ phu 

N-*+oo 



z = 



_ gtuo 

C 2 



l- e -p*> 2sinh(^)' 



c.q.f.d. 



La generalisation au cas a iV modes de bosons est immediate. Si on part du vide de 
l'espace de Fock, les etats coherents sout definis par 



e a |0 ) 



ou a est un nombre quantique qui decrit les fonctions sur la base desquelles l'espace 
de Fock est construit. Si on designe par z le vecteur de composantes z«, on a : 



a a z = z a z . 



La fonction de partition s'ecrit : 



Z ^ / Vz* a Vz a exp(-S[z*,z}), 

J z a (f3)=z a (0) 
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avec 



S[z*,z]= [ 
Jo 



13 

dr 



d 

Yl z *^ t )q^ z ^ t ) + H Zp(t),z*(t), z a (r), z p (t), z 7 (r), . . . ) 

a 

ou H est dans l'ordre normal. 
Approximations : 

• Perturbations (si H — H + V) 

• Phase stationnaire (utile, mais mal controle). 

4.3 Integrate de chemin et theorie de perturbations 

Supposons que H = H + V , ou H est quadratique dans les operateurs de bosons, 
et V est d'ordre superieur. La fonction de partition s'ecrit : 

JV 



Z= lim Y\T\ dZ * akdZak e- s (^ z ^) 



tv^+oo J- J- J- J- 27ri 

k=l a 

avec 

N ( (3 (3 

S = I ( Z °* Z ak - Z *k-lZ*ak) + J^H ({Z* ak , Z ak }) + —V ({z* ak , Z ak }) 

k=l \ a 

Dans une expression de ce type, V est simplement une fonction de variables com- 
plexes, et on peut done ecrire : 



- S ({ z ak> z <*k}) — p- S 0({Zak> z ak})p 



^ E %V({z* ak ,z ak }) 

e 



Dans cette expression, on peut developper l'exponentielle de V 

n=0 ' \fe=l 

On arrive done a une expression qui, pour iV donne, est un developpement en puis- 

sances de V. Par ailleurs, si H est quadratique, e~ 5 ° est de la forme e ij , 
ou on a regroupe les deux indices (a, k) dans un meme indice. Un terme d'ordre n 
du developpement se ramene done au calcul d'une integrate du type 



/ 



Vz* Vz Ziiz i2 . . . z in z* n . . . z* n e ^ , 
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ou les indices ii, . . . ,i n ,jn, ■ ■ ■ , ji sont ceux qui apparaissent dans la somme sur k 
elevee a la puissance n. 

Or, ce calcul peut etre effectue tres simplement a l'aide du theoreme de Wick qui 
stipule que : 



J Vz* Vz zaz i2 . . . z in z* n . . . z* j± e *<* 
f Vz* Vze w 



Demonstration : D'apres le resultat de l'integration sur les fonctions de variables 
complexes, on a : 

n<i* n Vz*Vze - ' Y-'^h^h 
G(J*, J) = ^-rjj = e^ 

J Vz* Vze « 



Si on derive la premiere expression par rapport a J* ± J* 2 . . . Jj n . . . Ji, il vient 

~Y z * M ij z j 



d 2n G 



fVz'VzznZe 



*jn ■ ■ ■ 



J=J*=0 



dJ* x . . . d.J n 

Si on derive la deuxieme expression, il vient 



-Y z t M H z j 

J Vz* Vze « 



d 2n G 



8 J*,... dJ n 



gr, 



J=J*=0 



81* 81* 

UJ il ■ ■ ■ UJ in 



x e< 



Y ■ J i M Tj lj 3 



j=j*=o 



Dans le calcul de cette derivee, les termes obtenus en derivant l'exponentielle sont 
nuls en J = J* = 0. Par ailleurs, pour obtenir une contribution non nulle, il faut 
que 



/ Vz* Vz z a z i2 . . . z in z* n . . . z* n e « 



~Y z i M ij z j 

J Vz* Vze « 



{k h l = l,...,n} = {i h I = 1, . . . ,n} 

-T,z*M, jZj 

— = E M ^--- M Av 



c.q.f.d. 



Le calcul des termes de la serie de perturbation repose sur l'utilisation systematique 
de cette identite sur la version discretisee des integrates de chemin. Dans ce contexte, 
le theoreme de Wick est une simple identite algebrique. 
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Remarque: On peut aussi faire une approximation de phase stationnaire. Si la solu- 
tion qui minimise Taction correspond a z = z* = 0, le developpement asymptotique 
autour du col est equivalent au developpement perturbatif. Par contre, si le min- 
imum est a z, z* 7^ 0, l'approximation de phase stationnaire conduit a un autre 
developpement, dit en nombre de boucles, et peut se reveler une approximation 
nettement meilleure. 



4.4 Integrale de chemin et fermions 

La representation de la fonction de partition et des fonctions de Green en termes 
d'integrale de chemin peut etre etendue au cas des fermions, mais il s'agit d'une 
extension non triviale. Ceci est du au fait que les operateurs creation et annihilation 
anticommutent. Du coup, la construction des etats coherents n'est pas immediate. 

Considerons en effet un etat coherent de bosons \z) et deux operateurs destruction 
a a et ap. On a : 



D'apres les regies de commutation des operateurs creation et annihilation de bosons, 
ceci doit donner 0, ce qui est compatible avec le fait que z a zp = zpz a pour z a , zp G C. 

Supposons qu'on reussisse a construire des etats coherents de fermions tels que 
a a \z) = z a \z), ap\z) = zp\z). D'apres les regies d'anticommutation, on doit avoir 



Les coefficients z a , zp ne peuvent etre des nombres complexes. lis doivent appartenir 
a une algebre de Grassman. 

Algebre de Grassman 

Une algebre de Grassman est definie par la donnee de n generateurs ip a , a — 1, . . . , n. 

• Ces generateurs anticommutent : ip a ipp + ippip a = 0. 

• En particulier, i\) 2 a = 0. 

• L'algebre comprend les combinaisons lineaires a coefficients complexes de tous 
les polynomes de degre au plus egal a 1 dans une variable donnee. 





(a a ap + apa a )\ z) = 

Z a Zp = —ZpZ a . 



70 CHAPITRE 4. INTEGRALE DE CHEMIN ET PROBLEME A N-CORPS 



Si le nombre de generateurs n = 2p est pair, on peut definir la conjugaison par 
une bijection d'un sous-ensemble de p generateurs vers le sous-ensemble des p 
restants. Le conjugue d'un generateur ip a est note ip a . 



Si a est un nombre complexe, (aip) = a*ip. 



• ^1^2 = 1p2*Pl- 

• On definit de fagon formelle l'integration par : 

J dip = J di> = 
dip ip — 1 / dip tp — 



• Pour un polynome de degre > 1, il faut anticommuter les variables pour 
ramener chaque variable ip a au contact de dip a . 

Exemple : J d^i d^ 2 V^V^i = J d^i V'sV'i = ~^3- 



Etats coherents 

Pour construire un etat coherent de fermions, on associe a chaque a a un generateur 
ip a et a a+ le conjugue tp a . Uespace de Fock generalise contient les combinaisons 
lineaires d'etats de l'espace de Fock avec des coefficients dans l'algebre de Grassman 
associee a ces generateurs. 

On impose que les variables de Grassman et les operateurs a, a + anticommutent, et 
on definit la conjugaison du produit par 

Par analogie avec le cas des bosons, on definit 

\ip) = e a |0 ). 



Proposition 4.1. 



i^)=n( i -^)i°: 



Demonstration : e ^ aCLa = 1 — ip a a^ puisque les termes d'ordre superieur, qui con- 
tiennent ip^, sont nuls. Par ailleurs, ip a at commute avec ippap. 
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Proposition 4.2. a a \ip) = ifj a \ijj) . 

Demonstration : 

car a a commute avec iftpap. Mais 

a Q (l - ip a a+) |0) = a Q |0 ) +ip a a a a+ |0) 

=o 

= jp a (l - t* a o+) |0 ) puisque ^ = 0. 

=^ a a \ip) = ip a \i/j) 

c.q.f.d. 

De meme (-^1 = (0|e « = (0|e«^ a " est un bra propre de a+ : (i>\a+ = 
Produit scalaire 

Considerons deux etats coherents construits a partir de deux algebres g et g' de 
generateurs ip a et ip' a . 

W> = (0|IJ(l + ^«o«)(l-^)|0> 

= (0| 11(1 + ^X^)10) 

a 

= (0|J](l + ^(l-a+a a ))|0) 

a 



Relation de fermeture 



| nd^d^e"?^|^)(^| = l. 



Demonstration: Appelons A cet operateur. On se propose de demontrer que 
(ai... a„|A|/3i . . . /? m ) = ( ai . . . a n \pi . . . f3 m ) 
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ou on a introduit les bras et kets propres de l'espace de Fock dermis par: 

\(3i . . . (3 m ) = . . . OgjO), {oti...a n \ = (0\a an . . . a ai . 

Les produits scalaires de ces etats avec les etats coherents sont donnes par: 

(ai . ..a n \ip) = (0\a an ...a ai \ip)= ...ij> ai 

et 

{■ip\(3i . . . (3 m ) =i/>p 1 ...i/>i3 m - 

Ainsi, 

(«! . . .a n |A|/?i ...0 m )= / JJd^ad^a JJ(1 - i>ai>a)lpa n • ■ -Ipai^pi ■ ■ ■ ^/3 m • 

Pour un etat 7 donne, les integrates qui peuvent intervenir sont 



J dip-y dip 7 (1 — ipjipy) 



( ] 




fil 


) ^ 7 1 








hi 




I t j 







Pour que le resultat soit non nul, il faut que chaque etat 7 soit vide ou occupe 
simultanement pour les cci . . . a n et /3i . . . (3 m . 

(a 1 ...a n \A\p 1 ...p m ) = (a 1 ...a n \p 1 ...p m ). 
Integrale gaussienne 



= / J_J_d^d^e u 1 

J i=i 

= det H e l - J 



Demonstration 



+00 I L~d Yi ±J -i]Y] 



m=0 



Dans l'integrale, seuls les termes contenant ipiipi . . . ipn^n donnent une contribution 
non nulle. Du coup, seul le terme m = n donne une contribution non nulle. Par 
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ailleurs, chaque terme non mil a comme coefficient un terme de det H. De plus, 
chaque terme apparait n\ fois puisqu'on peut choisir l'un des facteurs parmi n 
facteurs, le suivant parmi (n — 1) facteurs, et ainsi de suite. Enfin, on peut s'assurer 
que chaque terme est affecte du signe qu'il aurait dans le determinant 

1(0, 0) = det H. 

Exemple : 

J d^dtpe-^ = J d^d^(l - = a. 

c.q.f.d. 

Le calcul de I(fj, rf) se fait par la methode habituelle apres le changement de variables 
4>i = A - H^rjj. 



Trace d'un operateur 



( «i . . . On\lp ) = 1^ an . . . lp ai et ( IplP! . . . P m ) = Ipfo . . . 1pp m 

(a 1 ...a n \ij)(ij\p 1 ...p m ) = (ij\p 1 ...p m )(a 1 ...a n \ij)(-ir n 

car les variables de Grassman anticommutent. Si m et n ont la meme parite, 
(-1)" = (-l) m , et on a: 

(a 1 ...a n |V)(V'|A---/5m> = (-l) m (^\Pi...[3 m }(a 1 ...a n \^) 

= ( -iplPi-.-Pm )(a 1 ... a n \ip). 

Du coup, si A est un operateur qui conserve la parite du nombre de particules, on 
a : 



Tr A = ^2(n\A\ 



n 



I 



n 
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Fonction de partition 

Cette formule permet d'ecrire la fonction de partition sous la forme : 
Z= f nd^d^e"?^ Q (-^|e-^|^). 

J a 

En utilisant la formule de Trotter, cela se reecrit : 



n ™^ bosons 

1 la,r _ 2-ki 

n ajT d^ a (r)d^ a (r) fermions 

1 bosons 
-1 fermions 



Avec les notations 

et f = 

on peut ecrire 

f ^7^, -/o' dT fe^a(T)^:^a(T)+H(^ a ^ a ) 

Z — Vlp a V%p a e V a 

ou les ^ a sont des variables complexes pour les bosons et des variables de Grassman 
pour les fermions. 

Comme pour les bosons, cette formulation conduit de fagon tres economique a la 
theorie de perturbations. Par contre, si on essaie de faire une approximation de 
phase stationnaire, le champ correspondant au minimum est en general en dehors 
de l'espace de Fock physique car il s'exprime dans l'espace de Fock generalise avec 
des composantes dans l'algebre de Grassman. Pour pouvoir faire des approxima- 
tions de phase stationnaire significatives pour les fermions, il faut utiliser une autre 
representation de la fonction de partition en termes d'integrale de chemin a l'aide 
de champs auxiliaires. Pour cette nouvelle forme d'integrale de chemin, la represen- 
tation a l'aide des variables de Grassman est encore un outil tres utile, meme s'il 
n'est pas indispensable. 



Transformation de Hubbard- Stratonovitch 

Considerons un hamiltonien d'electrons en interaction decrit par 

a,P a/375 
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oil a,(3, . . . incluent le site et le spin. 

^Wc+c+c^ = - ^ V a/3lS c+ c+c 7 c s 

af55 a/3-yS 
=>• H — ^2 K a pf>al3 + ~ ^2 Votj3-y5Pa~iPl38 



2 

a,/3 af}^8 



avec 



K a 8 = t a 5 — ~ ^2 et p aj = C^C^ 







Le point de depart de la transformation de Hubbard-Stratanovitch est 1'identite 
gaussienne 



Si on fait le changement de variables Xi = ^ AijXj, on a : 

j 

l xAx = xA^x 

{[dxie « " = (2 7 r) iV/ Vdet Ae 

i 

La fonction de partition du hamiltonien precedent s'ecrit : 

p p 



Appliquons l'identite gaussienne au dernier terme avec 

A' 1 = -V mn , Xi = <7i, ji = p m , 771= (a, P), p al = 

C'est legitime puisque les variables p m commutent. Ainsi, 

2 12 PmVmnpn , C -| — r d(Xrri 2 12 cr mVmn&n 12 ^rnVmnPn 



2 12 PmVmnpn r- — f -| — r d(T m 2 12 

e m - n = vdet V / — ==e m - n 

m v 



La fonction de partition s'ecrit : 

D 8 - ( i A 

a/37* p "70 \ /3<5 / 



Z= Vip a Vtp a Vae a ^ so 
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ou le prefacteur 



Vdet V 



est inclus dans la definition de la mesure Va. 



On peut maintenant faire l'integration sur les variables de Grassman : 



Z = Va det 



d_ 



1 06 J 




C'est une nouvelle forme d'integrale de chemin. II faut faire les remarques suivantes : 

1. Les cr(r) sont des fonctions qui ne sont meme pas continues car il n'y a pas de 
terme cinetique. 

2. On peut faire une approximation de phase stationnaire. Le resultat a un 
sens physique. Comme on peut neanmoins faire plusieurs decouplages de 
l'hamiltonien en un terme quadratique et un terme biquadratique, on peut 
obtenir differentes approximations. 

Exemple cfcjCiCk = —5jkcfci + cfckcjci (voir plus haut) peut aussi etre de- 
compose comme : 



On peut demontrer que le premier decouplage correspond a l'approximation 
de Hartree, le deuxieme a l'approximation de Hartree-Fock. 

3. Cette representation et un certain nombre de variantes, notamment en termes 
de variables auxiliaires de type Ising, sont a la base de simulations numeriques 
de type Monte Carlo Quantique. 



ce qui conduit a l'hamiltonien 




avec 




Chapitre 5 



Integrate de chemin et theorie quantique 
des champs 



L'application de la methode des integrates de chemin a la theorie quantique des 
champs est une extension directe de la reformulation de la mecanique quantique par 
Feynman. Etant donne un lagrangien decrivant un champ, la quantification dite 
canonique consiste a identifier des champs conjugues et a postuler des relations de 
commutation entre ces champs conjugues. La methode des integrates de chemin 
consiste a ecrire les fonctions de correlation du champ comme les derivees d'une 
fonctionnelle generatrice, et d'ecrire cette fonctionnelle generatrice sous la forme 
d'une integrate de chemin sur tous les champs classiques a l'aide de Taction classique 
sous la forme 



L'objet de ce chapitre est de donner une introduction sommaire a cette methode en 
explicitant trois aspects : 



• Inequivalence de la quantification canonique et de la quantification par l'integrale 
de chemin dans le cas du champ de Klein-Gordon. 

• Inequivalence avec la formulation en termes d'integrale de chemin sur les etats 
coherents dans la limite non relativiste. 

• L'utilite, pour ne pas dire la necessite, de cette formulation dans le cas des 
theories de jauge. 



Z 
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5.1 Champ de Klein-Gordon 



La generalisation relativiste de l'equation de Schrodinger d'une particule libre sans 
spin proposee par Klein et Gordon consiste a partir de E 2 = p 2 c 2 + m 2 c 4 



iH d \ ^ ( ft \ ^ 

111 J <p(f f ,t)=\jVj V (f,t)+m 2 c 2 V (f,t), 



c dt / 

soit, en prenant k — c — 1, et en definissant le d'Alembertien par 



(□+m 2 Mx)=0 (x = (f,t)). 



C'est 1' equation de Klein-Gordon. 

Son interpretation en termes de fonction d'onde a une particule est problematique: 

1. C'est une equation du second ordre. Du coup, pour connaitre le mouvement 
ulterieur, il ne suffit pas de se donner la valeur de la fonction d'onde a l'instant 
initial, mais il faut aussi se donner la valeur de sa derivee. Une particule 
relativiste aurait done plus de degres de liberte qu'une particule non relativiste. 
On peut montrer que cela correspond a la charge de la particule. 

2. Comme l'energie n'est fixee que par la relation E 2 = p 2 c 2 + m 2 c i , il y a des so- 
lutions d'energies positives et negatives donnees par: E = ±\Jp 2 c 2 + m 2 c A . 
Les solutions d'energie negative doivent etre interpretees comme des anti- 
particules. 

3. L'equivalent de la densite en mecanique quantique non relativiste, 

ne satisfait pas d'equation de conservation. La quantite 

ih f dip dip*\ 

qui satisfait l'equation de continuite 

dp 



at ■ v ->' = ° 



avec 

h 



n'est pas toujours positive: ep(f, t) doit etre interprets comme la densite de 
charge. 
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En mecanique quantique relativiste, on doit l'interpreter en supposant que tp(x) est 
un champ. Cette equation est une equation de propagation pour un champ classique. 
Pour quantifier cette theorie, on peut proceder de deux fagons. 



Quantification canonique 

Cette equation peut etre derivee d'un principe de moindre action. Considerons en 
effet Taction : 



d 4 xC [(p(x),d^(x)} 



^ = («9 ,V) 
d» = (d , -V) 



avec 



L est la densite lagrangienne. 
Cette action est stationnaire si 



dip 
~dt 



- (Vip) 2 - m 2 ip 2 



Mais 



6<p(x ) 
5S[<p] d 2 ip 



5ip(x ) 

Le champ conjugue de (p est defini par 

U(x) 



dt 2 



+ Atp — m 2 (p. 



dC 



c.q.f.d. 



d(d (p) dt 

La densite hamiltonienne est la transformee de Legendre de C : 
soit 

H= 1 -n 2 + 1 -(^) 2 + 1 -m 2 v 2 . 

La theorie est quantifiee en supposant que les champs cp et II sont des operateurs 
qui satisfont les regies de commutation : 

' [if(ft),U(f',t)} = \5{f-f') 
[<p(rt),<p(r',t)}=0 

{ [n(ft),n(f',t)] = o. 
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On passe en transformee de Fourier : 

, - <1 : V. ; 

<p(r,t) 



avec uj^ = Vk 2 + m 2 . 



= 5 3 (k-k') 



H = somme d'oscillateurs harmoniques. 

Les fonctions de correlation du champ ip(r : t) s'obtiennent a partir de cette decom- 
position de Fourier. 



Quantification par l'integrale de chemin 

Le postulat de cette fagon de quantifier est le suivant : les fonctions de correlation 
du champ dans le vide sont donnees par la fonctionnelle generatrice 

Noter l'identite formelle avec la fonction de partition de Landau- Ginsburg. 
Plus precisement, 

et des formules equivalentes pour les fonctions a points. 

On peut demontrer explicitement pour le champ libre que cette expression est equiv- 
alente a celle obtenue par la quantification canonique. 

Etapes du calcul : 

1. Integration par partie pour les derivees temporelle et spatiales — > integrale 
gaussienne en ip. 

2. Calcul de l'integrale gaussienne. 

3. Calcul de la derivee seconde. 
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5.2 Limite non relativiste 



Pour prendre la limite non relativiste, on fait un changement de variables, et on 
definit : 



Jmt 



^(f,t) + ^n(f,t) 



■0 + (r, t) = e 
On verifie aisement que 



imt 



[0(f,t),0 + (f',t)] = <5 (3) (r-f') 
^(r,0,^(r',0] = [0 + (r,t),0+(f',t)] = 0. 



Cette transformation s'inverse en 



ip(r, t) = -h= [ip(r, t)e-' imt + ijj + (f, t)e imt ] 



n(f,t) = -i^f [ip(f,t)e 



-imt 



0+(f,t)e im *] . 



La densite lagrangienne 



s'ecrit : 



C = Ud ip - n 



<"5 



, , 90 90 + , 



^(V0+)(V0) + e 2imt [ ] + e- 2imt [ }. 
2m 



La limite non relativiste s'obtient en supposant que ip(r,t) et -0 + (f, t) varient lente- 
ment par rapport a e ±2imt . En effet, pour une particule non relativiste, la fonction 
d'onde <p(r,t) doit varier comme e~ tmt que multiplie une fonction qui varie lente- 
ment avec le temps. Ainsi, if)(r,t) et -0 + (r, t) contiennent une composante qui varie 
lentement avec le temps. Du coup, on peut negliger ces termes dans l'integrale de 
chemin. Apres une integration par parties, il vient : 



/ 



VipVip* exp 



J dt J dr0*i^ — 0" 



dt 



2m 



■0 + termes de source 



Cette fois, on a une derivee premiere par rapport a t. Cette expression est la 
meme que celle qu'on obtiendrait pour des bosons libres a l'aide d'operateurs champ. 
Inversement, on pourrait retrouver l'integrale de chemin en termes d'etats coherents 
en passant en transformee de Fourier et en utilisant une formulation generalisee avec 
une integrate sur VII. 
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5.3 Integrate de chemin et theories de jauge 



De meme que la quantification par l'integrale de chemin est beaucoup plus directe 
pour obtenir les fonctions de correlation que la quantification canonique, c'est une 
methode tres puissante dans le contexte des theories de jauge. Ceci dit, la presence 
d'un degre de liberte de jauge cree une difficulty supplemental. En quantification 
canonique, il y a des moments conjugues qui s'annulent. Considerons par exemple 
le lagrangien du champ electromagnetique 

C = - 1 -F^F, u -fA II , 



avec 



pvjx 



d 



(./" (./"../'). ./•„ (./". -./,). 0" = ^-...)- 

II ne contient pas la derivee par rapport au temps de A = ip (potentiel scalaire) 
> II" 0. 

Si Ton continue dans cette voie, A commute avec toutes les composantes Ilj =>- 
c'est un nombre. On a done perdu la covariance manifeste de la theorie. 

Dans la quantification par l'integrale de chemin, on pourrait etre tente d'ecrire 
simplement la fonctionnelle generatrice 



J VA^ exp i J d 4 xC {A^)-fA^ 



ou C (A^) peut s'ecrire, apres une integration par partie, 



□ 



9nv 



A v . 



Cette integrale a l'air d'une integrale gaussienne. Malheureusement, l'operateur 
" ~^\dnv — d^dv n'a pas d'inverse. II donne si on l'applique a un champ du type 
c^A, ou A est une fonction quelconque de x. Pour pouvoir utiliser la quantification 
par l'integrale de chemin, il faut fixer la jauge, ou plus exactement integrer sur 
les classes d'equivalence de configurations de A. En pratique, cela conduit a un 
lagrangien modifie pour le champ electromagnetique 



Co 



-A* 
2 



□ 



A v 



(5.1) 



ou jCs est le lagrangien de Stiickelberg, a est un parametre et d^A^ est un terme 
fixant la jauge. On peut verifier que les quantites physiques sont independantes de a, 
et que les resultats sont les memes que ceux obtenus par la quantification canonique. 
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Dans le cas de la chromodynamique quantique, ou Ton est confronte a des theories 
de jauge non abeliennes, cela impose d'introduire des champs fermioniques fictifs 
appeles "fantomes de Fadeev-Popov". 



Chapitre 6 

Appendices mathematiques 



A Integrates gaussiennes 

1. Z(j) = J dxe~^ Ax2+lx , Z(0) = J dxe~^ Ax 
Calcul de Z(0) : 

Z(0) 2 = J dx J dye 



Calcul de Z(j) 



= 2tt 

2tt 
~A 



r+oo 




/ rdr e" 


'0 




r+oo 




1 due 


'0 




- e -Au- 


+oo 






-A _ 






iA(x 2 +j/ 2 ) 



Am 



2(0) = J T 



1 , 2 ■ 1 , 

■-Ac + jx = --A 

= -ix 

2 



Z{j)= J dx'e-^ Ax ' 2 e^ 2/A = e&Z(0). 
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2. N variables : Soit A une matrice symetrique definie positive, et considerons 
l'expression 



/ 



N 

Z(j) = / Y[dxi exp 
J i=i 



j AT TV 

— — xiA^xj + 



r=l 



V 




On effectue le changement de variables x = x'+A 1 j. Comme x = i'+ j A" 
on en deduit aisement: 

1* 1* 1* 

— xAx +' jx = — x' Ax' H — j A 1 j 



2 2 

Z(j) = Z(0)e*' 



avec Z(0) = y J^dxj exp ^— - xAx^j . 



A est symetrique done diagonalisable. Soit R la transformation orthogonale 
(R l R = t) qui la diagonalise 



A =* i?L>i?, D 



di 

\ 

d N 



d i >0Wi( y A> 0). 



i? orthogonale det i? = 1. 

Si on fait le changement de variables x' = Rx, on a done 



N /i* \ 

j^J dxj exp I — xAx J 

i=l ^ ' 



/ n ^ ex p ( _ 2 x> 

N 

/n^ 

7 i=i 



exp ( -^dix'i 



-1/2 



m n/2 n* = < 2ir > 



7V/2_ 



Vdet A 



^■ ( 27r ) iV/2 i'M- 1 ," 



(0.1) 
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B Integrates gaussiennes et variables complexes 



1. 



dz dz* 
2vri 



-z'Hz+J'z+Jz* = H e R + 



7(0) = / 



d(Re z)d(lm z) n _ H ^ 



/rdr d6 
IT 

-L 



-nr- = I 2rdre- Hr2 



y = z- 



due~ Hu 
J 

H * 



-Hu 



-H 



z = y + 



o 

J_ 
H 



1 

H 



-z*Hz + J*z + Jz* = 



|J|2 |J|2 |J|2 

-y*Hy - y*J - J*y - i-±. + J* y + J y * + J-J- + J-J- 



Ul 



/(J) = ¥ e 



_lp|J| 2 /ff 



iV 

2 /n 



dz* dZj - z +Hz+J+z+z+J 



2vri 



hermitien positif, z et J vecteurs. 



hermitien =>- il est diagonalisable par une transformation unitaire R + R = 11 
telle que H = R + DR. 

Le jacobien de la transformation associee vaut 1. En effet, y = Rz. Posons 
y = yi+ k/2, z = zi + iz 2 , R = R\ + ii?2- 

R unitaire =>- RR + = 1 
(R 1 +iR 2 )( t R 1 -i t R 2 ) = 11 

i?i t R 1 + R 2 l R 2 = 1 
R 2 — R\ *i?2 = 
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Posons Y 



Z = 



Z\ 



, Z 2 

y = Rz <^> 



yi = Riz x - R 2 z 2 
y 2 = R 2 zi + R\z 2 



soit 



Y 



Demontrons que la matrice 



X 



R\ — i?2 

R 2 Ri 



Ri —R 2 
R 2 R± 



Z. 



est orthogonale. 

R2 Ri J \ — t R2 t R\ 

P t P 1 p (n P i P P (D 

r? it? p t p p t p 1 p t p 

11 

11 / 

c.q.f.d 

Si on fait le changement de variable z — > y dans l'integrale, le jacobien est egal 
a 1, et il vient : 

AT 



'(oh/ii^-^iU 

i=l i 



1(0) = (det H) 



-1 



Enfin, si J 7^ 0, on commence par le changement de variables 

y = z-H- 1 J z = y + H- 1 J 



Hz + J + z + z + J = 



-{y + + J + H^)H{y + if _1 J) + J+(y + H^J) + (y + + J + H^)J 
-y + Hy + J+H'iJ 




Remarque: Cette identite reste vrai meme si H n'est pas hermitien (demonstration 
non donnee), la seule condition etant la positivite de sa partie hermitienne pour 
assurer la convergence de l'integrale. 
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C Approximation de la phase stationnaire 



Considerons l'integrale 



dte iA/W 



et supposons que f'(t ) = 0. La contribution a cette integrale venant des points 
proches de t peut etre obtenue en developpant f(t) autour de t : 

/(*) = /(*o) + ilZ ^f"(t ) + 0((t - t ) 3 ). 
Si on ne retient que les termes d'ordre 2, on trouve 

I ~ e 



iA/(t„) /' +OO dte iA/"(t )(^fl^_ 
</— oo 



II s'agit d'une integrale gaussienne, et elle vaut 



2tt 



-iA/"(*o) 




C'est l'approximation de la phase stationnaire. Dans quelle mesure s'agit-il d'une 
bonne approximation ? On peut demontrer que l'integrale / consideree comme 
fonction de A possede un developpement en serie en puissances de ^, i.e. 



/(A) 



2i7r jWo) f 1 + ^1 + ^ + , 



A/"(*o) 



A A 2 



Par consequent, le terme qu'on a calcule domine dans la limite A — > +oo. Mal- 
heureusement, cette serie n'est pas convergente. elle est simplement asymptotique. 
Posons en effet 



/(A) 



2ivr 



3 iA/(to) 



A/"(*o) 

et ecrivons le developpement de F(X) sous la forme 



F(\) 



Alors 



n=0 



N 



n=0 
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Par contre, le rayon de convergence de la serie (en j) est egal a 0. 

La propriete fondamentale des developpements asymptotiques est que les premiers 
termes donnent une approximation tres precise si A est petit. 

Le meme type d'idee peut etre utilise pour une integrate du type 

»+oo 



r 

/(A) = / 

J — c 



dte" A ^. 



La contribution dominante vient du point ou / est minimum 

/'(*o) = 
f"(t ) > 




Dans ce cas on parle de methode du col (ou de Laplace). Les proprietes de conver- 
gence sont plus faciles a etablir dans ce cas. 



D Elements de calcul fonctionnel 



Une fonctionnelle est une fonction de l'espace des fonctions vers IR (ou C) 



F : {/} - IR 

/ -> F[f]. 



Autrement dit, une fonctionnelle associe un nombre a une fonction donnee. 

•t 2 



rt2 

Exemple: S[x] = / dt L(x, x,t) est un fonctionnelle. 
Jti 



Le concept de derivee peut etre etendu aux fonctionnelles, mais il s'agit d'une 
generalisation non triviale. 

Derivee de G&teaux : la derivee de Gateaux de F dans la direction v est une 
fonctionnelle definie par : 

ou v est une fonction. 
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Exemples : 

. F[f] = f( Xl )...f(x n ) 

D v F[f] = lim [/(^i) + ev ( x i)] ■ ■ ■ IfM + ev(x n )] - f(x 1 ) . . . f(x n ) 

Pour prendre la limite, il suffit de garder les termes lineaires en e 

D v F[f] = v(x 1 )f(x 2 )...f(x n ) + f(x 1 )v(x 2 )...f(x n ) + 
h f(x- I) ...f{x n - 1 )v(x n ). 

. F[f] = V(f(x)) 

e-+0 e 

Mais 

V(f(x) + ev(x)) = V(f(x))+ev(x)V'(f(x)) + 0(e 2 ) 
=► D v F[f] = V(f(x))v(x). 

. F[f] = JdxG(f(x)) 

D v F[f]= J dxG'(f(x))v(x). 

. F[f] = V(f'(x)) 

D v F[f] = V\f{x))v'{x). 

. F[f] = JdxG(f'(x)) 

D v F[f] = J dxG'(f(x))v'(x). 

D'apres la definition de D v F[f\, il est clair que l'application qui a v associe D v F[f] 
est une forme lineaire sur l'espace des fonctions, done une distribution 



D v F[f] = J dxv(x 



SF[f] 

sf(xy 



La derivee fonctionnelle de F par rapport a / est la distribution definie par cette 
equation. En pratique, s'obtient a partir de D v F[f] en remplagant v(x) par 

oj(y) 

5(x-y). 
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Exemples : 

. F[f] = V(f{x)) 

On verifie bien stir que 

A TP\ f 1 C 

d y v (y)T7T\ = / dyv(y)V'(f(x))6(x-y) 



Sf(y) 

= v(x)V'(f(x)). 



• F[f] = J dxG(f(x)) => ^ = G"(/(y)) 
. F[/] = V(f(x)) 

Attention ! 

| dy ^(/'(^MyWs-i/) = ^'(/'(x))t;'(x) 

mais 

| d^/^M*-?/) = -V ff (/'(y))/ ff (y) 
car 5' est impaire. 

. F[f} = J dxG(f'(x)) 
SF[f\ 



Sf(y) 



dxG'(f(x))S'(x-y) = -G"(f(y))f"(y). 



D.l Proprietes 



La derivee fonctionnelle ainsi definie possede les proprietes habituelles de la derivee 



Sf(x) 1 6f(x) 6f(x) 
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D.2 Derivees d'ordre superieur 

D V2 D Vl F[f] est bien definie puisque D Vl F[f] est une fonctionnelle. C'est manifeste- 
ment une forme bilineaire en v\ et t> 2 , et elle s'ecrit de fagon generale 

5 2 F[f] 



D V2 D Vl F[f] = J dxiwi(xi) J dx 2 v 2 (x 2 ) 



Sf(x 1 )Sf(x 2 ) 



5 2 F\f] 

ou T77 — tttt — 7 es t un produit de distributions par rapport aux variables x± : x 2 , . . . 
of(xi)5f(x 2 ) 

On definit de meme les ordres superieurs. 



D.3 Developpement de Taylor 



F[f + v] = F[f} + J dxivM 



1 /"j / X /"j , s S F if] 

+ 2j dxMXl) J dX ^hf(x 1 )5f(x 2 ) + - 
Demonstration : on considere la fonction 

git) = F[f + t n ] 

, (t) = ^FUHt + eH-Flf + t,] 

= D v F[f + trj}. 
Or, si g {n) {t) = D„ . . . D r , F[f + trj], il est evident que 

N v ' 

n fois 

g^)(t) = D 11 ...D v F[f + tr j }. 

S v ' 

n+1 fois 

On en deduit que 

g(t) = g(0) + 1 D v F[f] + t.Dr,D v F[f\ + ... 
Pour t = 1, il vient : 

F[f + V] = F[f] + D v F[f] + ^D v D v F[f] + ... 
Par definition de la derivee fonctionnelle, 

S (n) F[f] 



p v . . .D^ F[f] = J dxi 77(^1) . . . y da;„ 77(x„) 



fois 



Sf(xi) ■ ..6f(x n ) 



94 



CHAPITRE 6. APPENDICES MATHEMATIQ UES 



E Transformee de Fourier 



Soit Aij une matrice reelle telle que = A(i — j) ne depende que de i — j. On 
suppose qu'on a des conditions qux limites periodiques. Alors A^ est diagonalisee 
par une transformee de Fourier. Plus precisement, supposons que i — 1, . . . , N, et 
considerons la matrice unitaire 



Soit 



1 - 27rm 

U mn = -^e 1 f n , m,n = 1, . . . , N. 



N 



A m n ~ (U AU 1 ) mn 

= UmiAjjUjn 



h3 



- J_ V Jqi a . Jet' j J 9 - -W m 

1,3 



N ^ lJ I q ' = f n 

l-Y^e^-^A^'^ 

i,3 



NS qq' cette somme est independante de i 
(A invariante par translation et con- 
ditions aux limites periodiques). 



Definitions : La transformee de Fourier de la fonction A(i — j) est definie par 

A(q) = J2 A (r) e iqr . 



La matrice A = (UAU x ) s'appelle la transformee de Fourier de la matrice A. Si A 
est periodique, A est diagonale et est donnee par: 



A mn — 8 mn A ( q — ~j^ m ) • 
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La transformee de Fourier inverse s'ecrit : 

E\r x ~A u ■ 
u im ^hnn u nj 



1 ^— \ ■ 2TTim ~ ( 2lT 

= — > e n A [ q — —m | e* * 



2irmj 



l^e-^U(9). 



Si on passe a la limite continue, on se limite a un intervalle de longueur 2n, 
par exemple [— n, n] 



Si A est tridimensionnelle, on integre sur la premiere zone de Brillouin (ZB) 
dont le volume est (27r) 3 



A(r t - r 3 ) = J^y 3 J zB dqe-^-^A (q) . 



Enfin, A- 1 = (Ay 1 . 
En effet, 



A- 1 = UA- l U- 1 

= UAU^ 1 = A 



=S> A- 1 = (A) \ 

L'ingredient essentiel est que la matrice de passage est la meme pour toutes 
les matrices invariantes par translation. 



• Consequence : 

- £E (*•>)" 



1 ^e^^ 



